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Abstract. This paper presents an algorithm to compute the value of 
the inverse Laplace transforms of rational functions with poles on ar- 
rangements of hyperplanes. As an application, we present an efficient 
computation of the partition function for classical root systems. 



1. Introduction 

The ultimate goal of this work is to present an algorithm for a fast com- 
putation of the partition function of classical root systems. We achieve this 
goal in somewhat more general terms, namely we develop algorithms to 
compute the volume of a polytope and its discrete analog, the number of 
integer points in the polytope. These formulas, in turn, are inverse Laplace 
transforms of certain rational functions, and our work can be viewed in these 
general terms. 

Let U be a finite-dimensional real vector space of dimension r. Denote 
its dual vector space U* by V. Consider a set of elements 

A = {ai,a 2 , • • • ,&n} 

of non-zero vectors of V. We assume that the convex cone C(A) generated 
by non-negative linear combinations of the elements on is an acute convex 
cone in V with non-empty interior. 

The elements £ in V produce linear functions u i— ► £(u) on the complexified 
vector space Uq- In particular, to the set A we associate the arrangement 
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of hyper planes 



N 



H C (A) : 



\J{u£U c \ tti(u) = 0} 



i=i 



in Uc and its complement 



C/ C (^l) := 



ueU c \ n«i(n) / 



i=i 



We denote by 7£_4 the ring of rational functions on Uc(A) with poles along 
TCc(A). Then each element (j) £ Ha can be written as P/Q where P is a 
polynomial function on r complex variables and Q is a product of elements, 
not necessarily distinct, of A. 

Our first aim is to present an algorithm to compute the value of the inverse 
Laplace transform of functions in TZ_\ at a point h £ V. In other words, we 
study the value at a point h € V of convolutions of a number of Heaviside 
distributions <fi i— > J °° (f)(toti)dt. The first theoretical ingredient is the notion 
of Jeffrey-Kirwan residues [14]. Going a step further, DeConcini-Procesi [12] 
proved that one can compute Jeffrey-Kirwan residues using maximal nested 
sets (in short MNS), a combinatorial tool related to no-broken-circuit bases 
of the set of vectors A. 

The applications in view are volume computation for polytopes, enumer- 
ation of integral points in polytopes and, more generally, discrete or con- 
tinuous integration of polynomial functions over polytopes. Indeed, Szenes- 
Vergne [21], refining a formula of Brion-Vergne [7], stated formulae for the 
volume and number of integral points in polytopes involving Jeffrey-Kirwan 
residues. 

Consider the polytope 



As a function of h, the volume of II^/i) is a piecewise-defined polynomial. 
The chambers of polynomiality in the parameter space V are polyhedral 
cones. 

Our programs are extremely efficient for computing the volume of the 
polytope II_4(/i) when A is a classical root system. An important fact is that 
our algorithm can work with formal parameters, thus giving the polynomial 
volume formula for Tl^h) when h runs over a particular chamber. 

For an analogous theory for integral-point enumeration, we have to as- 
sume that the a\ are vectors in a lattice V%. For h € V%, the function N^(h) 
which associates to the vector h the number of integral points in Hj^(h), 
that is the number of ways to represents the vector h as a sum of a certain 
number of vectors aj, is called the (vector )-partition function of A. For 
example for E-2, given a vector (hi, /12) with integral coordinates we would 
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like to compute the number NB 2 (h) of vectors (xj) G such that 

Xl ( ) +X2 ( 1 ) +X3 ( -1 ) +X4 ( 1 ) = ( h 2 

As a function of /i, the number N^(h) of integral points in H^,(h) is a 
piecewise-defined quasipolynomial, and again the chambers of quasipolyno- 
miality are polyhedra in V [19, 20]. 

In this paper, we describe an efficient algorithm for MNS computation for 
classical root systems. This algorithm for MNS gives rise to programs for 
Kostant partition function for the classical root systems A n , B n , C n , and 
D n . Again, our algorithm works with a formal parameter h that is assumed 
to be confined to a particular chamber. 

These calculations are valuable because partition functions play a fun- 
damental role also in representation theory of semisimple Lie algebras g. 
Indeed, partition functions arise naturally when we want to compute the 
multiplicity of a weight in a finite-dimensional representation or the tensor- 
product decomposition of two representations, both being basic problems 
to understand characters of representations. Cochet [9] has obtained very 
efficient algorithms for both these problems in the case of A n , implementing 
results of [2]. See also a forthcoming paper [10] for multiplicities computa- 
tion in all the classical Lie algebras using the results obtained in this paper. 
There is also a class of infinite-dimensional representations, the discrete- 
series representations, whose understanding is central for the general theory 
of admissible irreducible representations. The decomposition of such repre- 
sentations to a maximal compact subgroup of g is predicted by Blattner's 
formula, which is a partition function in which the roots involved are the 
so-called noncompact roots. 

We conclude by describing the way the paper is organized. Section 2 
introduces Laplace transforms and polytopes. In Section 3, we recall Jeffrey- 
Kirwan residues and its link with counting formulae. DeConcini-Procesi's 
maximal nested sets are described in Section 4, as well as how they are 
related to Jeffrey-Kirwan residues. Section 5 describes our general algorithm 
for MNS computations. Details of particular cases of the algorithm for the 
root systems A n , B n , C n and D n are examined in Sections 7-10. Finally 
comparative tests of our programs with existing softwares are performed in 
Section 11. 

A number of theoretical results on the function Nj^{h) when A is a subset 
of the system A n can be found in Baldoni-Vergne [1] (as, for example, the 
computation of the volume of the Chan-Robbins polytope). 

Computer programs for volume computation/integral-point enumeration 
in polytopes have only been implemented in the very recent past, most 
notably LattE [16, 17] and barvinok [6], both of which are implemen- 
tations of Barvinok's algorithm [3]. To the best of our knowledge, these 
two are the only general programs for volume computation/integral-point 
enumeration in polytopes. More specialized programs include algorithms 
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of Baldoni-DeLoera-Vergne for flow polytopes [2] and Beck-Pixton for the 
Birkhoff polytope [4]. 

Our programs have been especially designed for classical root systems, 
are faster than all actual existing softwares and can compute new examples 
that were not reachable by previous algorithms. Note in particular that our 
programs can perform computations for N^(h) for A n at least up to n = 10 
(11 coordinates vector). For B n , C n , T> n the algorithms are efficient at least 
up to n = 6. For our methods (as well as for LattE), the size of the vector h 
affects only little on the computation time. Recall that our methods can also 
calculate the multivariate quasi-polynomials h ^ N^(h) when h varies on 
a chamber, and as a particular case for a fixed h the function k ^ N^(kh) 
which is the Ehrhart quasipolynomial in k. 

2. Laplace transform and polytopes 

We start by briefly recalling the notations of the introduction, aiming to 
relate the inverse of the Laplace transform with various counting formulae 
for a polytope. A good introduction on this theme is the survey article [23]. 

2.1. Laplace transform. Let U be a finite-dimensional real vector space 
of dimension r with dual space V. We fix the choice of a Lebesgue measure 
dh on V. Consider a set 

A = {a>i,a 2 , ■ ■ ■ ,a n } 

of non-zero vectors of V. We assume that the set of vectors a« spans V. For 
any subset S of V, we denote by C(S) the convex cone generated by non- 
negative linear combinations of elements of S. We assume that the convex 
cone C(A) is acute in V with non-empty interior. 

Let V s ing(A) be the union of the boundaries of the cones C(S), where S 
ranges over all the subsets of A. The complement of V s i ng (A) in C(A) is by 
definition the open set C reg (A) of regular elements. A connected component 
c of C reg {A) is called a chamber of C(A). Figures 1 and 2 represent slices 
of the cones C(A^) and C(B^), where the dots represent the intersection of 
a slice with a ray M>o hence showing the chambers. Note that the cham- 
bers for B, r and C r are the same (as roots in B r and C r are proportional). 
In dimension 3, the root system A3 is isomorphic to D3. See [2] for the 
computation of chambers. Very little is known about the total number of 
chambers. On the other hand, given a vector h, it is easy to compute the 
equations of the chamber containing h. This was done in [2, 11]. We have 
incorporated this small part of the corresponding program in our programs 
for classical root systems. 

Table 3 represents the only numbers of chambers that have been computed 
(and the computation time). 

Consider now a cone C(S) spanned by a subset S of A and let p be a 
function on C(S). We assume that p is the restriction to C(S) of a polynomial 
function on V. By superposing such functions p, we obtain a space CV(V, A) 
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Figure 1. The 7 chambers for A3 
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Figure 2. The 23 chambers for 
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Figure 3. Number of chambers and computation time 



of locally polynomial functions on C(A). For / S CV(V,A), the restriction 
of / to any chamber c of C(A) is given by a polynomial function. 
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The Laplace transform L(f) of such a function / is defined as follows. 
Consider the dual cone C(A)* C U of C(A) defined by: 

C(A)* = {u G U | (h,u)>0 for all h G C(.4)} . 

Then for u in the interior of the cone C(A)*, the integral 

L(/)(u) = / e-^f(h)dh 
Jc(A) 

is convergent. It is easy to see that the function L{f) is the restriction to 
C(A)* of a function in 7?._4. (Recall that 1Zj± is the ring of rational functions 
P/Q on U where P is a polynomial function on U and Q is a product of 
elements of A.) It is easy [7] to characterize the functions L(f) on U arising 
this way. 

Let v be a subset of {1,2,..., n}. We will say that v is generating (re- 
spectively 6osic) if the set {oil \i £ u} generates (respectively is a basis of) 
the vector space V. 

Every basic subset is of cardinality r and we write Bases („4) for the set 
of basic subsets. Given a 6 Bases(„4), the associated basic fraction is 

In a system of coordinates (depending on a) on U where cti(u) = Ui (for 
i £ a), such a basic fraction is simply of the form 

1 

U\U2 ■ ■ ■ U r 

Define Q(U, A) C TZa as the linear span of functions ^ — — n-, where v is 

generating and n» are positive integers. The following proposition gives the 
characterization we were speaking of and is easy to prove: 

Proposition 2.1. [7] If f is a locally polynomial function on C(A), the 
Laplace transform L(f) of f is the restriction to C(A)* of a function in 
Q(U,A). Reciprocally, for any generating set v and every set of positive 
integers n{ > 0, there exists a locally polynomial function f on V such that 

TT 1 ~TV^ = I e- {k ' U) f( h )dh 

for any u in the interior of C(A)* . 

We define the inverse Laplace transform L _1 : Q(U,A) — > CV(V,A) as 
follows. For <f> G G(U, A), the function L~ l (f> is the unique locally polynomial 
function that satisfies 

<f>(u)= f e -< fc>tt >(rV)(/»)(tt 

Jc(A) 

for any u £ C(A)* . 

In the next sections, we will explain the relation between Laplace trans- 
forms and the enumeration of integral points of families of polytopes. We 
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will see in Section 4 that one can write efficient formulae for the inversion of 
Laplace transforms in terms of residues, whose algorithmic implementation 
is working in a quite impressive way, at least for low dimension. 

2.2. Volume and number of integral points of a polytope. In this 
subsection we consider a sequence 

A + = [Oil, C62, ■ ■ ■ , CCjv] 

of non-zero elements of A. We assume that each element a £ A occurs in 
the sequence; in particular N > n and the set A + spans V. 

Remark 2.2. In all our examples, the sequence A + will not have multiplic- 
ities, so that we will freely identify A + and A. 

We introduce now the notion of a partition polytope. 
We consider the space R N with its standard basis u>i and Lebesgue mea- 
sure dx. 

If x = Yli=i x i u i ^ ^ N w ith Xi > (1 < i < N) then we will simply write 
x > 0. 

Consider the surjective map A : R N — > V defined by A(ui) = a { and 
denote by K its kernel. Then K is a vector space of dimension d = N — r 
equipped with the quotient Lebesgue measure dx/dh. 

If h G V, we define 

U A +(h) = {x e R N \ Ax = h;x > 0} . 

The set H^+(h) is a convex polytope. It is the intersection of the non- 
negative quadrant in M. N with an affine translate of the vector space K. 
This polytope consists of all non-negative solutions of the system of r linear 
equations 

N 

^XiCti = h. 

i=l 

Remark 2.3. It might be appropriate to recall that any full dimensional 
convex polytope P in a vector space E of dimension d, defined by a system 
of N linear inequations 

P = {yeE\(u i ,y) + \ i >0} 

(where Ui G E* and Aj are real numbers), can be canonically realized as a 
partition polytope IT_4+ (h). Here A + is a sequence of iV elements in a vector 
space of dimension r = N — d. Indeed, consider the diagram 

E -U R N -A, v = R N /i{E) 

where i : y i-> J2iLi { u iiU) UJ i an< i A is the projection map R N — ► V. Let 
ai be the images of the canonical basis u>i of M. N . Define A + = [a\ , . . . , ajv] 
and consider the point h := A(%2iLi \^>i)- Then the polytope U^+(h) is 
isomorphic to P. Indeed, the points in H^+(h) are exactly the points x; t 
such that YliLi( x i ~ = with xi > 0. By definition of the space 
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V = R N /i(E), there exists y G E such that xi — Xi = (ui,y). As Xi > 0, 
this means exactly that (ui, y) + Xi > 0, so that the point y is in P. 

More concretely, to determine the partition polytope Ax = b starting 
from a polytope P given by Qy T > A (where Q is a N x d matrix whose i th 
row is given by a vector n« G P* and A G E^) we choose among the elements 
Ui a basis of E* . Thus after relabeling the indices and doing an appropriate 
translation, we may assume the inequations of the polytope P are given in 
the form 



C is a r x d matrix, 
A G E r , 

and y = (yi, . . . 6 



Vi > o 
2/2 > 

where 

Vd > 

Cy T + X T > 

Then the polytope P is isomorphic to the polytope defined by 

jx > 0| ^lx T = A T } 
where A is the r x N matrix given by 

-C I r 



A 



I r being the identity matrix, 
be the polytope defined by the system of in- 



Example 2.4. Let P C 

equalities: 

-xi + 1 > 0, 

-x 2 + 2 > 0, 

-xi - x 2 + 2 > 0, 

2x 1 +x 2 -l > 0. 

Choosing the basis u\ = (1,0), u 2 = (0, 1) and using the translation y\ 
—x\ + 1 and i/2 = —x 2 + 2 we can rewrite the system as: 

f yi>o 



V2 > 

y\ 

1J2 



c 



+ 



-1 



where C 



> 



1 1 

-2 -1 



Therefore P is isomorphic to 

n^ + (M = {y=(2/i,...,2/4)GE 4 ,t/>0 

with h 



-yi -V2 + V3 = 

2yi + 2/2 + 2/4 = 



-l 

3 



We continue with our review. If h is in the interior of the cone C(A), then 
the polytope (/i) is of dimension d. It lies in a translate of the vector 
space K, and this translated space is provided with the quotient measure 
dx/dh. 
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Definition 2.5. We write vol^+(h) for the volume of Il^+(h) computed 
with respect to this measure. 

Suppose further that V is provided with a lattice Vz and that 

A + := [ai,a 2 , ■ ■ ■ ,a N ] 

is a sequence of non-zero elements of Vz spanning Vz, that is, Vz = YliLi ^ a i- 
In this case, the lattice Vz determines a measure dzh on V so that the 
fundamental domain of the lattice Vz is of measure 1 for dzh. However, for 
reasons which will be clear later on, we keep our initial measure dh. We 
introduce the normalized volume. 

Definition 2.6. The normalized volume volz,^+ (h) is the volume of II 4+ (h) 
computed with respect to the measure dx/dzh. 

Remark 2.7. The reason for keeping our initial dh is that the root systems 
B r ,C r , D r live on the same standard vector space V = W, where the most 
natural measure is the standard one. This measure is twice the measure 
given by the root lattice in the case of C r and D r . 

If vo\(V/Vz,dh) is the volume of a fundamental domain of Vz for dh, 
clearly volz j\+(h) = vol(V/Vz, dh)vol^+ (h) . 

Let now h G Vz- A discrete analogue of the normalized volume of II 4+ (h) 
is the number of integral points inside this polytope. 

Definition 2.8. Let N^+(h) be the number of integral points in U^+(h), 
that is the number of solutions x = (x±, . . . , xn) of the equation J2iLi x i a i = 
h where x\ are non-negative integers. The function h 1— > N_^+(h) is called 
the partition function of A + . 

We will see after stating Theorem 3.3 that the functions h 1— > vol^+(h) 
and h 1— ► N^+(h) are respectively polynomial and quasipolynomial on each 
chamber of C(A). 

The following formulae (see, for example, [23]) compute the Laplace trans- 
form of the locally polynomial function vol_4+(/i) and the discrete Laplace 
transform of the quasipolynomial function N^+(h). 

Proposition 2.9. Let u G C(A)* ■ Then: 

(1) / e-^Kol A+ (h)dh= 1 
JC(A) Ui=l a i( u ) 

(2) y e-^jy (/ i ) = 1 

heV z nC(A) l^il 1 > 

3. Jeffrey-Kirwan residue 

The aim of this section is to explain some theoretical results due to Jef- 
frey and Kirwan which are fundamental for our work. They described an 
efficient scheme for computing the inverse Laplace transforms in the context 
of hyper plane arrangements. 
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Let's go back to the space of rational functions 1Zj\.- It is Z-graded by 
degree. Of great importance for our exposition will be certain functions in 
7Za of degree — r. Every function in 1Zj± of degree — r may be decomposed 
into a sum of basic fractions f a (see Equation (1)) and degenerate fractions; 
degenerate fractions are those for which the linear forms in the denominator 
do not span V . Given a £ Bases(„4), we write C{a) for the cone generated by 
OLi (i £ a) and by vol(<r) > for the volume of the parallelotope X^i=i[P> M a i 
computed for the measure dh. Observe that vol(cr) = | det(cr)|, where a is 
the matrix which columns are the ctj's. Now having fixed a chamber c, we 
define a functional JK c (0) on TZ^ called the Jeffrey-Kirwan residue (or JK 
residue) as follows. Let 



(2) JK C (/ CT 




if C C C(cr) 
if cnC(a) 



By setting the value of the JK residue of a degenerate fraction or that of 
a rational function of pure degree different from — r equal to zero, we have 
defined the JK residue on TZ^- 

We may go further and extend the definition to the space TZ_a which is 
the space consisting of functions P/Q where Q is a product of powers of 
the linear forms aj and P = J2h=o ^k is a formal power series. Indeed 
suppose that P/Q 6 TZ_a where we may assume that Q is of degree q, and 
P = J2T=o Pk is a formal power series with P^ of degree k. Then we just 
define 

JK C (P/Q) = JK c (P q - r /Q) 
as the JK residue of the component of degree — r of P/Q. In particular if 
4> G TZa an d h G V, the function 



k 



k=0 



is in TZ_a and we may compute its JK residue. Observe that the JK residue 
depends on the measure dh. 

Let's now make a short digression that should clarify why JK residues 
compute inverse Laplace transforms. For u 6 C(A)* we have: 



-L- f e-^dh = f a (u). 



In other words the inverse Laplace transform of f a computed at the point 
h € C(A) is vo \( a ) Xa(h), where \u is the characteristic function of the cone 
C(a). We state this as a formula: 

^x.(») = £-'(/.)(*). 

Since the JK residue can be written in terms of basic fractions, the fol- 
lowing theorem [14] is not surprising: 
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Theorem 3.1 ( Jeffrey-Kirwan) . If <f> € TZa> then for any h 6 c we have: 

{L- 1 4>){h)=m c (e< h ' •><£). 

Assume that ^ : U — > f/ is a holomorphic transformation defined on a 
neighborhood of in £7 and invertible. We also assume that aj(F(u)) = 
otj(u)fj(u), where fj(u) is holomorphic in a neighborhood of and fj(0) / 0. 

If is a function in TZ^, the function ^*<p(u) = <p(^{u)) is again in TZ^. 
Let Jac(^) be the Jacobian of the map The function Jac( x I / ) is calculated 
as follows: write ^(u) = U2, ■ ■ ■ , u r ), . . . , ty r (ui, U2, ■ ■ ■ , u r )). Then 

Jac(\P)(u) = det((^-*j)ij). We assume that Jac(*)(u) does not vanish 

at u = 0. For any in 7£_4 the following change of variable formula [1], 
Theorem 45, which will be useful in our calculations later on, holds: 

Proposition 3.2. The Jeffrey- Kirwan residue obeys the rule of change of 
variables: 

JK C (^) = JK c (Jac(^)(**0)). 

We conclude this section by recalling the formula for N^+(h). 

Consider the dual lattice Uz = {u G U\ (u, Vz) C Z} and the torus 
T = U /Uz- Choosing a basis {u±, . . . , u r } of Uz we may identify T with the 
subset of U defined by the fundamental domain for translation by Uz- 



with < tj < 1. 



Every element g in T = U /Uz produces a function on Vz by h i— > 
e (h,2nV^iG) ^ w ]j ere we denote by G a representative of g € U/Uz- 1 For 
u G Bases(„4) we denote by T(cr) the subset of T defined by 



T(a) = 5 G T 



; (a,2^^TG> = x for alla G a j _ 



This is a finite subset of T. In particular if er is a Z-basis of Vz, then T(a) 
is reduced to the identity. More generally, consider the lattice TLo generated 
by the elements a in a. If p is an integer such that 7Lo C pVz, then all 
elements of T(a) are of order p. 

For g £ T and h £ Vz, consider the Kostant function F(g, h) on U defined 

by 

e (/i,2?r v^-TG+u) 

(3) F(g, h){u) = —f, — -. 

For example when g = 0, 

p {h,u) 

F(0,h)(u) 



n£i(i-e-<°"«>)' 



1- We prefer to denote the complex number i by ^/^T because we use i for many indices. 



12 M.W. BALDONI, M. BECK, C. COCHET, M. VERGNE 

The function F(g, h)(u) is an element of Indeed if we write 

1(g) = j j | 1 < { < N, e -<«i,27rv^TG) = j j ^ 

then 

p (h,u) 

(4) /*)(«) = e <MW=TG> 

Hie/(9) \ ai ' u / 

where ip 9 (u) is the holomorphic function of u (in a neighborhood of zero) 
defined by 

1 h9( u \= TT x TT - 

If c is a chamber of C(A), the Jeffrey-Kirwan residue JK c (F(g, h)) is well 
defined. 

The following theorem is due to Szenes-Vergne [21]. If the set A is uni- 
modular (that is, each a £ Bases (A) is a Z-basis of Vz), it is a reformulation 
of Khovanskii-Pukhlikhov Riemann-Roch calculus on simple polytopes [15]. 
For a general set A, this refines the formula of Brion-Vergne [7] . 

Theorem 3.3. Let c be a chamber of the cone C(A) andc its closure. Then: 
(1) For h £ c we have 

( e^>-> \ 

vol z A + (h) = vol(V/Vz, dh) JK C 



nN I - 

t=l a i J 

(2) Assume that F is a finite subset ofT such that for any a G Bases(„4), 
we have T(a) C F. Then for h G Vz l~l c, we /iaue 

(/») = vo\(V/V z , dh) JK C (F( 5 , fc)). 

Observe that the right-hand side of (2) does not depend on the measure 
dh, as it should be. 

Let us explain the behavior of these functions on a chamber c. By defi- 
nition, a quasipolynomial function on a lattice L is a linear combination of 
products of polynomial functions and of periodic functions (functions con- 
stants on cosets h + pL where p is an integer). We now show that the 
normalized volume vol z _4+(/i) is given by a polynomial formula, when h 
varies in a chamber c, while N A + (h) is given by a quasipolynomial formula 
when h varies in V% l~l c. 

The residue vanishes except on degree — r, so that 



.imi <<*,«>; vn£i <<*,«>; 

and as expected the normalized volume is a polynomial homogeneous func- 
tion of h of degree iV — r on each chamber. 
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Now if A is unimodular then the above defined set F is {0}. Hence the 
number of integral points in the polytope II 4+ (h) satisfies 



N A+ (h) = vol (V/Vz, dh) JK C 



,{h,u) 



nf =1 (i- e -(^)). 



and is a polynomial of degree N—r whose homogeneous component of degree 
N — r is the normalized volume. More precisely, write 

nili(i - e-<«*.«>) " nli <<*,«> x nL(i - e - <Qi,M> ) 

where 

i=i <a».«> 



nii(l-e-<«"«>) 



fc=0 

is a holomorphic function of 11 in a neighborhood of with tpo(u) = 1. 
Consequently 

iY4 + (/i) = vol(y/y z ,d/i)JK c -x^^fe(n) 



A;=0 



,u) 



The unimodular case applies to the root system A r . 

Finally in the non-unimodular case (for example for the root systems 
B r , C r , D r ) the set F is no longer reduced to {0}. Let us denote by 
ip 9 (u) = ^2k^o^k( u ) the series development of the holomorphic function 
V> 9 appearing in formula (4). Then we see that 3K c (F(g, h)) equals 

/ , Jh,u) \ 

(6) JK C e (h,2,V^iG)^_ r{u) 



h d^)|-»-*)' 'I Hie/to) K«>^ (,j ' 



If g is of order p, the function /i e (h,2ir^/~^G) j g cons t an t on each coset 

h + pVz of the lattice »Vz, while the function /i JK c ( iii^ , v ib 9 (u) ) 

is a polynomial function of h of degree (/(fl 1 )! — r — k. Thus the function 

(7) N M (h) = vol (V/Vz, dh) JK c (F(g, h)) 

g&F 

is given by a quasipolynomial formula when h varies in the closure of a 
chamber. Note that its highest degree component is polynomial and is the 
normalized volume as expected. 
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Example 3.4. Let us compute the normalized volume and number of in- 
tegral points for the root system B2, that is for A + = 62 = {ei,e2,ei + 
62; ei -62}- Fix a chamber c and an integral vector h = (hi, /12) in the cone 
0(82)- Observe that the root lattice is Ze\ © Ze2 and vol (V/Vz, dh) = 1 for 
the measure dh = dh\dh2- Then the normalized volume equals 

J_ JK ( (h 1 u 1 + h 2 u 2 ) 



2! \uiU2(ui + u 2 )(ui — 112) 
Note that 

«? 1 

+ 



Ullt2(«l + U 2 )(U1 — t/ 2 ) U2(ltl+U2) (Ul + U 2 )(U1 — U2) 

(and similar quotients of u\ and 1/1U2), so that the normalized volume is 

1 TT ^ / hj h\ + 2h 1 h2 + h 2 2 h\ 

■ JK C — + 



2 \U2{UI+U 2 ) (U\ + U 2 ){ui - l/ 2 ) Ul(t/l+lt2) 

There are three chambers, namely Ci = C({e2, ei + 62}), C2 = C({ei, ei+e2}), 
C3 = C({e\ — e2, ei}) (see Figure 4). Now let us compute the Jeffrey-Kirwan 
residues on the chambers. As 



JK 
JK 



ci \u 2 {u 1 +u 2 ) ) J Kc 2 ( (ui+u 2 )(ui-u 2 ) ) 2' 

1 I 1 TTy I 1 \ _ 1 



c 2 yui(tll— U2)/ ' C3 \ (ui+«2)(Ml— « 2 ) / 2' 

we obtain 

vo\{U B2 (h)) = -hj if/iGCi, 

voi(n S2 (/*)) = 1(^ + ^)2-1^ jf/je^ 

voi(n B2 (/*)) = \{hi + h 2 f if // e c 3 . 



Note that the formulae agree on walls Ci PI C2 and C2 fl C3. 

For the number of integral points, we first note that F = {(0, 0), (1/2, 1/2)}. 
Consequently iVg 2 (h) is equal to the Jeffrey-Kirwan residue of f\ = F((0, 0) , h) 
plus / 2 = F((l/2, 1/2), h). We rewrite the series fj (j = 1, 2) as fj = 
ft x e ulhl+U2h2 /uiU2(ui + n 2 )(t*i - u 2 ) where 

7*1 V>2 7*1 +U2 U\ — U2 



1 — e -«i l — e~" 2 1 — e~( ni+ " 2 ) 1 — e~( Ul ~" 2 ) ' 
1 + e - " 1 1 + e~ u 2 1 _ e -(ui+u 2 ) i_ e -(«i-« 2 ) v ' 

Using the series expansions x * e - x = 1 + \x + y^x 2 + 0(a; 3 ) and l+ x e - x = 
\x + 0(x 2 ), we obtain that the number of integral points is the JK residue 
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e, + e 9 




of 



Figure 4. The 3 chambers for B 2 



i*i(l + f/ii + \h\) i | + /11/12 + |/i 2 + 5 /ti , u 2 {\hl + \h 2 ) 



+ 



+ 



« 2 (U1 — U 2 )(ui + U2) («1 — U2){ll\ + U 2 ) U\{U\ - U 2 )(ui + U 2 ) 

1 

+ ( — 1 )' ll+/l2 i 



(ui + n 2 )(ui - u 2 ) 
{l + lh x + \h{) l + 2(hi + h 2 ) + \h 2 2 + h 1 h 2 + \hl (U| + U2 



u 2 (m + u 2 
+(— 1 )^ 1+?i2 



+ 



(mi - u 2 )(ui + u 2 ) 



U\(u\ + u 2 ) 



We then obtain: 

N B2 (h) = 

N B2 (h) = 
N B2 (h) = 



(ui + u 2 )(ui - u 2 ) 



1 + 2*1 + 2*? 



1 



if /t € Ci, 
1 



1 



1 



7 



\hi+h2 



if h € c 2 , 

Jfc? + \h x h 2 + ±h 2 2 + h 1 + h 2 + 7 - + (-l) 11 ^ 1 - 
if h e C3. 



1 



Note that the functions Nb 2 agree on walls, and the formulae above are valid 
on the closures of the chambers. 

Our general method to implement Theorem 3.3 for root systems is more 
systematic and will be explained in the course of this article. 
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Remark 3.5. Combining (6) and (7), we can see that the quasipolynomial 
character of the integral-point counting functions Nj^ stems precisely from 
the root of unity in (6). Furthermore, we will see in Lemmas 8.2, 9.1, and 
10.1 that for root systems of type B, C, and D, these roots of unity are 
of order 2, as in the above example for £>2- (For root systems of type A, 
(5) shows that N% is always a polynomial.) Let us record the following 
immediate consequence: 

Corollary 3.6. The integral-point counting functions Ns r ,Nc r ,Nx> r are 
quasipolynomials with period 2. 

Remark 3.7. The partition functions N_^ r , Nq t , Nq t , Nx> r can be inter- 
preted as (weak) flow quasipolynomials on certain signed graphs [5]. The 
polynomiality of Nj^ r follows immediately from this interpretation and a 
unimodularity argument; the fact that the quasipolynomials Ng r , Nc r , Nj) r 
have period 2 follows from a half-integrality result of Lee [18]. 

Remark 3.8. In the case where A + is an arbitrary sequence of vectors in Vz, 
the straightforward implementation of Theorem 3.3 above is of exponential 
complexity. Indeed we make a summation on the set F, which can become 
arbitrarily large. Barvinok uses a signed cone decomposition to obtain an 
algorithm of polynomial complexity, when the number of elements of A + 
is fixed, to compute the number N^+(h); the LattE team implemented 
Barvinok's algorithm [16, 17] in the language C. Our work will be dealing 
either with volumes of polytopes, where the set F does not enter, or with 
partition function of classical root systems, where the set F is reasonably 
small. Then we obtain a fast algorithm, implemented for the moment in 
the formal calculation software MAPLE. This algorithm for these particular 
cases can reach examples not obtainable by the LattE program. 

In next Section 4 we will give the basic formula for JK C involving max- 
imal proper nested sets, as developed in [12], and iterated residues. These 
formulae are implemented in our algorithms. 



If / is a meromorphic function of one variable z with a pole of order less 
than or equal to h at z = then we can write f(z) = Q(z)/z h , where Q(z) 
is a holomorphic function near z = 0. If the Taylor series of Q is given 
by Q(z) = YlT=o ( IkZ k , then as usual the residue at z = of the function 
f( z ) = ^2T=o ( lkZ k ~ h is the coefficient of 1/z, that is, qu-i- We will denote 
it by res 2= o/(z). To compute this residue we can either expand Q into a 
power series and search for the coefficient of z^ 1 , or employ the formula 



We now introduce the notion of iterated residue on the space IZ^. 



4. A FORMULA FOR THE jEFFREY-KlRWAN RESIDUE 
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Let v = [ai, a 2 , • • • , a r ] be an ordered basis of V consisting of elements of 
A (here we have implicitly renumbered the elements of A in order that the 
elements of our basis are listed first). We choose a system of coordinates 
on U such that oii(u) = Uj. A function <ft 6 TZa is thus written as a ratio- 
nal fraction (j){ui,u 2 , ■ ■ ■ ,u r ) = '"'„ r j where the denominator Q is a 
product of linear forms. 

Definition 4.1. If € IZa, the iterated residue IreSj?(</>) of for z7 is the 

scalar 

lresp((f)) = res Ur=0 res Ur _ 1= o • • • res ul= o^(ui, u 2 , ...,u r ) 
where each residue is taken assuming that the variables with higher indices 
are considered constants. 

Keep in mind that at each step the residue operation augments the homo- 
geneous degree of a rational function by +1 (as for example res x= o(l/xy) = 
1/y) so that the iterated residue vanishes on homogeneous elements 4> £ TZa> 
if the homogeneous degree of (f> is different from — r. 

Observe that the value of Ires^((/>) depends on the order of V. For example, 
for / = l/(x(y - x)) we have Tes x=0 res y=0 (f) = and reSj, =0 res x=0 (/) = 1. 

Remark 4.2. Choose any basis 71, 72, . . . , 7 r of V such that ® 3 k=1 a,j = 
©k =1 7j for every 1 < j < r and such that 71 A72 A - • • A7 r = a.\ f\a 2 A - • -Aa r . 
Then, by induction, it is easy to see that for (f> G TZa 

res ar= o • • • res ai= o0 = res 7r=0 • • • res 71=0 </>. 
Thus given an ordered basis, we may modify a 2 by a 2 + ca>i, . . . , with 
the purpose of getting easier computations. 

The following lemma will be useful later on. 

Lemma 4.3. Let V = [01,0:2, • • • ,a r ] and fp = ^r 1 /3 . be a basic fraction. 
Then the iterated residue Ires^/g) is non zero if and only if there exists a 
permutation w of {1,2, ... ,r} such that: 

(3 w{2) € MaielRas, 

(3 w ( r ) £ Rai • • • ®Ra r . 

Definition 4.4. Let V = [a\, a 2 , . . . , a r \ and let Uj = ctj{u). Choose a 
sequence of real numbers: < e\ < e 2 < ■ ■ ■ < e r . Then define the torus 

(9) T(v) = {u £ Uc I \uj\ = ej, j = 1, . . . ,r}. 

The torus T(v) is identified via the basis ctj with the product of r circles 
oriented counterclockwise. The sequence [e\,e 2 , . . . ,e r ] is chosen so that 
elements a q not in (B J k=1 Raj do not vanish on the domain {u £ Uc \ \uk\ < 
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< k < j ; \ui\ = €i,i = j + l,...,r}. This is achieved by choosing 
the ratios ej/ej+i very small. The torus T(u) is contained in Uc{A) and 
the homology class [T(j?)] of this torus is independent of the choice of the 
sequence of the ordered ej [22]. 

Choose an ordered basis e\, e2, • • • , e r of V of volume 1 with respect to the 
measure dh. For z G Uc, define Zj = (z, ej) and dz = dz\ A dzi A • • • A dz r . 
Denote by det(j?) the determinant of the basis a±, ct2, ■ ■ ■ , a r with respect to 
the basis e±, e2, • • • , e r . 

Lemma 4.5. For 4> £ T^A> we have 

deTM reSa ^° " " " reSai = 00 = v^hy /t ( .) 

Thus, as for the usual residue, the iterated residue can be expressed as 
an integral. 

We now introduce the notion of maximal proper nested set, MPNS in 
short. 

De Concini-Procesi [12] prove that the set of MPNS is in bijection with 
the so-called no broken circuits bases of A (with respect to a order to be 
specified). This is helpful as the JK residue can be computed in terms of 
iterated residues with respect to these bases. 

If S is a subset of A, we denote by (S) the vector space spanned by S. 
More generally if M = {Si} is a set of subsets of A, we denote by (M) the 
vector space spanned by all elements of the sets Si. We say that a subset S 
of A is complete if S = (S) H A or in other words if any linear combination 
of elements of S belongs to S. A complete subset S is called reducible if we 
can find a decomposition V = V\ © V2 such that S = Si U S2 with Si C V\ 
and S2 C V2. Otherwise S is said to be irreducible. 

Definition 4.6. Let X be the set of irreducible subsets of A. A set M = 
{7i, I2, ■ ■ ■ , Ik} of irreducible subsets of A is called nested if, given any sub- 
family {Ii, . . . ,I m } of M such that there exists no i, j with Ii C Ij, then 
the set Ii U ■ ■ ■ U I m is complete and the elements Ij are the irreducible 
components of Ii U I2 U • • • U I m . 

Example 4.7. Let E be an r + 1-dimensional vector space with basis 
{i = 1, . . . , r). We consider the set 

K-r+i = {ei — ej \ l < i < j < r + I}. 

These are the positive roots for the system A r . The irreducible subsets 
of IC r+ i are indexed by subsets S of {1,2, ... ,r + 1}, the corresponding 
irreducible subset being {e^ — ej\i,j G S, i < j}. For instance the set 
S = {1, 2, 4} parametrizes the set of roots given by {e\ — e2, e2 — e^, e\ — e^}. 

A nested set is represented by a collection M = {Si, S2, ■ ■ ■ , Sk} of subsets 
of{l,2,...,r + l} such that if Si, Sj G M then either Si n S*j is empty, or 
one of them is contained in another. 
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Definition 4.8. A maximal nested set (in short MNS) M is a nested set 
such that for every irreducible set 1 of A the set MU{I} is no longer nested. 

A maximal nested set has exactly r elements [12]. 

Assume now that A is irreducible, otherwise just take the irreducible 
components. Then every maximal nested set M contains A. Let I\, I2, ■ ■ ■ , 
Ik be the maximal elements of the set M\A. We see that the vector space 
spanned by (Ji) © (I2) © • • • © (Ik) is of codimension 1 [12, Proposition 1.3]. 

Definition 4.9. A hyperplane H in V is A- admissible if it is spanned by a 
set of vectors of A. 

Thus if M is a MNS, the vector space (M\A) is an admissible hyperplane 

H. 

Definition 4.10. Let A be irreducible and let H be a „4-admissible hyper- 
plane. All MNPS's such that (M \A) = H are said attached to H. 

Therefore to classify maximal nested sets (MNS) for an irreducible set A 
we proceed by running over the set of ^.-admissible hyperplanes, as described 
in Figure 5. 

• Take a hyperplane H spanned by a set of vectors of A. 

• Break A l~l H into irreducible subsets I\ U I2 U ■ ■ ■ U Ik ■ 

• For each irreducible U construct the set {M[, . . . , M£.} of maximal 
nested sets for ij. 

• Set Ci = {l,...,ki}. 

• A maximal nested set is then given by the union M\ x U M^ 2 U • • • U 

U {.4} where c\ G C\, . . . , c& € and all of them are obtained 
by letting Cj vary. 

Figure 5. Building of all MNSs attached to an ^l-admissible 
hyperplane H 

The whole algorithm will be described in detail in Figure 7, Section 5. 

We describe now the notion of maximal proper nested set of A. 

Fix a total order on the set A. For example, we can choose a linear 
functional ht on V so that the values ht(aj) are all distinct and positive. 
Thus the value ht(a) is larger if a is deeper in the interior of the cone. 

Let M = {S\, S2, ■ ■ ■ , <Sfc} be a set of subsets of A. In each Sj we choose 
the element ctj maximal for the order given by ht. This defines a map 9 
from M to A. 

Definition 4.11. A maximal nested set M is called proper if 9(M) is a 
basis of V. We denote by V(A) the set of maximal proper nested sets, in 
short MPNS. 

If M = {ii, I2, ■ ■ ■ , I r } is a maximal nested set, we associate to M the 
list [6(1^), . . . , 0(li r )\ using the total order on the elements 9(M); that is 



20 M.W. BALDONI, M. BECK, C. COCHET, M. VERGNE 

we have ht^ijj) < ht(0(I i2 )) < ••• < ht(0(J ir )). Observe that, if ^4 is 
irreducible, for every maximal nested set, h r is always equal to A and 9(Ii r ) 
is the highest element of A. We will often implicitly renumber our elements 
in M such that ht(0(/i)) < ht(0(I 2 )) < • • • < ht(0(I r )). 

So we have associated to every maximal proper nested set M an ordered 

basis 9(M) = [a±, a 2 , ■ ■ ■ , a r ] of elements of A. In all implementations, we 

> 

calculate 0(M) from a MNS M with the procedure ThetaMNS(M). We denote 
by vol(M) > the volume of the parallelepiped ^^ =1 [0, ljo;^ with respect 
to our measure, and by C(M) = Yll=i ^>o&i C C(A) the cone generated by 
9{M). 

If v is a regular element of V, let 

(10) V(v, A) = {M£ V{A) | v G C(M)}. 

The set V(v, A) depends only of the chamber c where v belongs. We are 
now ready to state the basic formula for our calculations. 

Theorem 4.12 (DeConcini-Procesi, [12]). Let c be a chamber and let v G c. 
Then, for (f> G TIa> we have 

We will use also the corresponding integration formula. 

Each maximal proper nested set M G V(v,A) determines an oriented 

cycle T (9{M)j contained in the open set Uc(A), as described in Defini- 
tion 4.4. 

Definition 4.13. Let c be a chamber. Define the oriented cycle: 

H(c) = Yl si S n ( det (W)) [T (Hm)) ■ 

MeP(v,A) 

The following integral version of Theorem 4.12 will be useful. 
Theorem 4.14. Let c be a chamber. Then for (j) G TZa we have 

JK c (0) = ~ )=- I ^{z)dz. 
(2vr v / -l) r Jh(c) 

The following example should help clarifying the notions introduced. 

Example 4.15. We consider the set /C4 of positive roots for ^3 (see Fig- 
ure 6) defined by 

lC 4 = {e i -e j \l<i<j< 4}. 
We let V be the vector space generated by the elements in /C4. Then V 
has dimension 3 and we write an element of V as 

a = aiei + a 2 e 2 + a 3 e 3 - (01 + a 2 + as)e4- 
We consider the height function defined by 

ht(ei - e 2 ) = 10, ht(e 2 - e 3 ) = 11, ht(e 3 - e 4 ) = 12. 
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<H + «3 > ai > 



Figure 6. Hyperplanes for ^3 with a = a\e\ + 02^2 + ^363 — 
(ai + a 2 + a^)ei 

This choice gives the following order on the roots: 

ei — e2 < e2 — e3 < e3 — e4 < ei — e3 < e2 — e4 < ei — e±. 

Take a hyperplane H in V spanned by two linearly independent elements 
of /C4. Therefore it is the kernel of a linear form ai, where In is a 

proper subset of {1, 2, 3, 4}. The set of complementary indices gives the same 
hyperplane. Thus each admissible hyperplane partitions the set of indices 
{1,2,3,4} in two sets Z\ and Z2, where Z\ := {i G Ih} and Z2 is the set 
of complementary indices. In our example we have 7 choices of admissible 
hyperplanes corresponding to the following partitions: 

^ = {[1,2,3], [4]}, ff 2 = {[1,2,4], [3]}, fZa = {[1,3,4], [2]}, 
H A = {[2,3,4], [1]}, ff 5 = {[1,2], [3,4]}, tf 6 = {[1,3], [2,4]}, 
F 7 = {[1,4], [2,3]}. 

Now observe that if the hyperplane Hi already contains the highest root 
ex — e4 then it cannot lead to a maximal proper nested set. Indeed we must 
get a basis if we add the highest root to a set of vectors contained in Hi. 
Thus H2, H3, H7 can be excluded. It remains to consider the hyperplanes 
Hi, H4, i/5, Hq. 

Hyperplanes Hi and H4 give rise to two MPNSs each, while H$ and Hq 
give rise to only one. So we obtain a list of 6 maximal nested sets (as 
described in Example 4.7, we identify an irreducible subset I with a subset 
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S of [1,2,3,4]): 

Mi = {[1,2], [1,2,3], [1,2,3,4]}, M 2 = {[2,3], [1,2,3], [1,2,3,4]}, 

M 3 = {[2, 3] , [2, 3, 4] , [1, 2, 3, 4]}, M 4 = { [3, 4] , [2, 3, 4] , [1, 2, 3, 4] } , 

M 5 = {[1, 3], [2, 4], [1, 2, 3, 4]}, M 6 = {[1, 2], [3, 4], [1, 2, 3, 4]}. 

5. Search for maximal proper nested sets adapted to a vector: 

the general case 

Given a vector v in the cone C(A), we describe how to search for all 
maximal proper nested sets belonging to V(v,A), without enumerating all 
MPNS. 

We use as height function a linear form that is positive and that takes 
different values on all elements ctj, and consider the total order it induces. 
Let H be an ^.-admissible hyperplane in V, that is, a hyperplane spanned 
by a set of vectors of A. Then the cone C(AnH) generated by the elements 
of A belonging to H is a cone with non-empty interior in H. 

We have already seen that to list all the MPNS, we have to first list all 
admissible hyperplanes H and then find the irreducible components J\, J 2 , 
. . . , J s of Af) H. Then we choose a MPNS M, t := {If} for J u and define 
M = Mi U M 2 U M 3 U • • • U M s U {^4}. 

As we have seen in Example 4.15 we can discard some of the hyperplanes 
a priori, because they cannot lead to a maximal proper nested set. The next 
lemma examines the general situation. Let 9 be the highest element in A 
and H a hyperplane of A. 

Lemma 5.1. There exists a maximal proper nested set M £ V(v,A) at- 
tached to H , if and only if 9 does not belong to H and if v belongs to the 
cone generated by 9 and An H. 

Proof. The condition is necessary. Indeed v must belong to the cone gener- 
ated by the elements 9 (If) and 9, and all the elements 9 (If) are in A n H. 
Reciprocally consider the projection v — j0jj9, where u is the equation of 
the hyperplane H. This can be written as v\ © V2 © ■ ■ ■ © v s , where each V{ 
is in the cone C(Ji). Let now Mj G V(vi, Ji) be a MPNS in Jj. The element 
Vi belongs to C(9(Mi)). We can write 

s 

i=\ if a Mi 

with tf > 0. Thus we see that the collection Mi U • • • U M s U A is a maximal 
proper nested set in V(v,A). Moreover in this way we list all elements of 
P(v,A). □ 

Our search for maximal proper nested sets in V(v , A) will then be pursued 
by constructing all possible admissible hyperplanes H for which v is in the 
convex hull of C(A H H) and 9. We denote by Hyp(v,A) the set of such 
^.-admissible hyperplanes. 
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The following easy lemma lists some obvious conditions for the set Hyp(-u, -4). 
Let uh £ U be the normal vector to an „4-admissible hyperplane, meaning 
that H := {h £ V \ (u H ,v) = 0}. 

Lemma 5.2. If H £ Hyp(v,A) then H satisfies the following conditions: 

(1) (u H ,e)^o. 

(2) (u H ,v) x (u H ,9) > 0. 

Thus if a hyperplane H satisfies the above conditions we define 

pio) H {v) = v- — 9. 

[UH, tf) 

Hence to decide if H £ Hyp(f , A) we simply have to test if proj H (v) is in the 
cone generated by A n H, which is done by standard methods. Our search 
for the hyperplanes H £ Hyp(f , A) will also be considerably sped up by the 
following remark. 

Proposition 5.3. Let H be an A- admissible hyperplane. Let u £ U be a 
linear form on V which is non negative on An H and on 9. If (u,v) < 0, 
then H is not in Hyp(v,„4). 

Proof. Indeed if v was in the cone generated by A n H and 9, the value of 
u would be non negative on d. □ 

The point of this remark is that in classical examples of root systems, an 
a priori description of the ^l-admissible hyperplanes is available, together 
with the defining equations of the cone C(AnH). This condition will allow 
us to disregard right away many „4- admissible hyperplanes. 

Let us summarize the scheme of the algorithm in Figure 7. Recall that 
we have as input a vector v, and as output the list of all MPNS's belonging 
to V(v,A). 

We will explain our algorithm in more details for each classical root system 
(see Sections 7-10). 

6. Trees and order of poles 

Let M be a maximal nested proper set for the system A := {a±, a2, ■ ■ ■ , a n }. 
In our algorithms, we will need to take an iterated residue with respect to a 
basis 9{M) of a function of the form <fi = ~, where P is a polynomial 

function on U. It is thus important to understand the order of the poles 
of the function obtained after performing a certain number of residues. We 
also prove that the iterated residue associated to M depends only on the 
tree associated to M. 

We associate to a maximal nested set M a tree T as follows. Let M = 
{I±, . . . , I r } be a maximal nested set. The vertices of T are the elements of 
M and the oriented edges are determined from the reverse order relation by 
inclusion: the ends of the tree are irreducible sets with just one element and 
if A is irreducible, the base is the set A. A subset N of M will be called 
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check if v € C(A) 

for each hyper plane H do 

check if v and 9 are on the same side of H 

if not, then skip this hyperplane 

define the projection proj H (v) of v on H along 9 

check if proj^(-u) belongs to C(A n -ff); if not then skip this hyperplane 
write A H .ff as the union of its irreducible components I\ U ■ ■ ■ U i& 
write v as v\ @ ■ ■ ■ ® Vk according to the previous decomposition 
for each Ij do 

compute all MPNS's for Vj and Ij 
collect all these MPNS's for Vj and Ij 
end of loop running across ij's 
collect all MPNS's for the hyperplane H 

end of loop running across Ws 

return the set of all MPNS's for all hyperplanes 

Figure 7. Algorithm for MPNS's computation (general case) 

saturated if it contains all elements above elements of N in the tree order. 
Thus if ./V contains an element S, it contains all the elements S' of M which 
are contained in S. 

Example 6.1. The two MNSs named Mi and M5 described in Example 4.15 
can be rewritten respectively as 

[1,2] [1,3] [2,4] 

[1,2,3] and 

[1,2,3,4] [1,2,3,4] 

Lemmas 7.3, 8.6, and 10.4 describe the decomposition of A PI H in irre- 
ducible nested sets and lead to the following result: 

Proposition 6.2. Let T be the tree associated to an irreducible classical 
root system. Then T is a connected tree for which every vertex is adjacent 
to at most two other vertices. 

Lemma 6.3. Let M = {I±, I2, ■ ■ ■ , I r } be a MPNS. Here we have numbered 
our irreducible sets such that 6(Ii) < #(^2) < ■ • ■ < 0{I r ). Let k be an integer 
smaller than or equal to r. Then the set {I±, I2, ■ ■ ■ , Ik} is saturated. 




Indeed if two sets /, J belongs to M and I C J, then 9(1) < 9 (J). 
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Proposition 6.4. Let M = [I±, I 2 , ■ ■ ■ , I r ] be a maximal nested proper fam- 
ily. Let [I[,I' 2 , . . . , I' r ] be a reordering of the sequence [Ji, I2, ■ • ■ ,Ir\- We 
assume that this reordering is compatible with the partial order given by 
inclusion: if Ij C I' k then j < k. Let 

v=[B{h),e{h),...,e{i r )\ 

and 

Then we have Ires,? = Ires p . 

Proof. We prove this proposition by induction on r. 

If A is irreducible, then necessarily L r = I' r = A and I' 2 , ■ ■ ■ , Z^-il 
is a reordering of the sequence [L\,L 2 , ■ ■ ■ ,I r -i]- Furthermore the families 
{Ii, I 2 , . . . , I r -i} and {I[,I 2 , . . . , i^-i) are maximal proper nested sets for 
Aq = U^z\lj. The set „4o spans a codimension 1 vector space in V. 

To prove that IreSj? = Ires^, , it suffices to test it on basic fractions f a . Let 
a = {(3i, (3 2 , . . . , (3 r } be a basic subset of A. By Lemma 4.3, if Ires^/o- 7^ 0, 
then the set a n (Aq) is of cardinality r — 1, and there exists an element of a, 
say p r , of the form c6 + £ where £ belongs to (Aq), c is a non-zero constant 
and 9 is the highest element of A. Let 

v = [9(i 1 ),e(i 2 ),...,e(i r - 1 )} 

and 

4 = [^),^(/ 2 ),...,0(/;_ 1 )]. 

Then we have 

Ires^/ CT = -Ires ul J an{Ao) 

and 

Ires^-X = -Ires^/ CTn( ^ o) . 

We conclude by induction. 

When A is not irreducible, we write A = U* =1 J a where J a are irreducibles. 
We have V = ©* =1 (J a ). Every basic subset a of A is the union of basic 
subsets for the irreducible sets J a . Define 

u a = [9(I^),9(I i a %...,9(J a )} 

where [I* 1 , I* 2 , . . . , J a ] is the subsequence of irreducible sets contained in J a 
extracted (with conserving order) from the sequence [I±, I 2 , . . . , I r \. Simi- 
larly let 

^ = [«(C)-»(C) «W1 

where [L' a n , L' a t2 , . . . , I' a ] is the subsequence of irreducible sets contained in J a 
extracted from the sequence [![, I 2 , . . . , I' r \. Then, as the calculation takes 



26 



M.W. BALDONI, M. BECK, C. COCHET, M. VERGNE 



place with respect to independent variables, we have 

s 

Ires /7 (/ CT ) = Yl(lre Sl ;J ar]{ j a) ), 

a=l 

s 

Ires-,(/ CT ) = n( Ire V'J-n(J a >)- 

a=l 

Each of the vector space ( J a ) is of dimension less than r, so that by induction 
hypothesis Ires^ = Ires^, . This concludes the proof. □ 

Let us now consider partial iterated residues. To a set v of elements of 
A, we associate the vector space 

H v := {u G U | (a, u) = for all a G v}. 

A linear function a £ A produces a linear function on H u by restriction. 
If v := [ai, CC2, ■ ■ ■ j «fc] is a sequence of elements of A, the partial iterated 
residue 

Ires^ := res afe= o • • • res ai=0 ^ 

associates to a rational function (f) in a rational function on H u of the 
form 

G 

Hi=l,...,n;ai^0 a i n * 

where G is a polynomial function on H u and a is the restriction of a to H u . 
Let M be a MPNS and consider the tree associated to M. Given a saturated 
subset S of M, we can define the iterated residue with respect to this satu- 
rated set: we choose any order S := [I±, I2, ■ ■ ■ , Ik] on S compatible with the 
inclusion relation and define Iress := IreSj? with V = [9(Ii),6(l2), ■ ■ ■ , 6(Ik)]- 
With the same proof as for Proposition 6.4, this partial residue depends 
only on the set S. We denote by H$ the intersection of the kernels of the 
elements a for a G S. It is also the intersection of the kernels of the elements 
9(Ik), as the set v is a basic sequence in S. 
Let 4> be a function in Rj± of the form 

<? = ™ • 

n i= i <*» 

Let M be a MPNS and J±, J2, ■ ■ ■ , J s be elements of M. We consider the 
saturated subset S of M consisting of the elements of the tree strictly above 
Ji, J2, . . . , J s . The iterated residue Iress^ is a function on Hs- Denote by 
u a the restriction of the function 9{J a ) to H$- 

Proposition 6.5. The pole of the linear function u a in the iterated residue 
lress<j> is of order less than or equal to \ J a \ — dim (J a ) + 1. 

See Figures 8 and 9 for an application of the proposition. 
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Proof. Choose a vector space E such that 

V = {J 1 )®---@{J S )@E. 
Let B = U* =1 J a and C = A \ B. Write C := {fa, fa, ■ ■ ■ , A?} and 

</> = PX(/>iX02X---X</> s X(3 

with = iri^ and Q = nj-bv 

For /?,- G C, we write ^ = f% =1 (5) + 7^ with fa- G ( J;) and 7,- G £7. The 
element jj is necessarily non zero, as the set B is complete. Thus we write 

1 _ 1 

and the iterated residue is by definition 

(q , 00 / y^s m" 

P X (0! • • • S ) X JJ - E (-1)"^^ 
j=l 7 J fe=0 V , 

Here, when taking the residue, the elements jj are considered as constants 
and this sum is finite. 

Consider the subset M a of elements of M contained in J a . This is a 
MPNS for the set J a . Let j£ be the saturated subset of M a consisting of 
all elements of M a different from J a . Then has dim ( J a ) — 1 elements. If 
g = jj — ° a „ a , the iterated residue IreSj+<7 is a Laurent polynomial in u a . 

Now Iresg^ is a sum of products of residues of the form IreSj+<7 a where 

g a = 1 r Pa — and P a is a polynomial. Thus we obtain a Laurent polynomial 

in u a , a = 1, . . . , s (with coefficients rational functions on the vector space 
E*). Now the homogeneous degree of g a is greater than or equal to — |J a |. 
The number of residues we are taking is equal to dim(J a ) — 1. So we 
obtain a function of u a of homogeneous degree greater than or equal to 
— I J a \ + dim ( J a ) — 1. This means that the pole in u a is of order less than or 
equal to \J a \ — dim (J a ) + 1. □ 

Let us consider the MPNS whose tree representation is given by Figure 8. 
The orders of the poles of its nodes are given in Figure 9. 

Remark 6.6. In our program for calculating iterated residues for root sys- 
tems of type A r , we will reorder roots according to the tree order: we take 
the residue first with respect to the elements 9(Ik) appearing at the end 
of the tree in arbitrary order, and we remove these variables. Then we 
take the variables appearing at the end of the tree when we have removed 
these irreducible sets. Here an irreducible set I is indexed by a subset S of 
{1, 2, . . . , r + 1}. A subset S of cardinality 2, for example [1, 3], corresponds 
to the irreducible set with one element (here e\ — e%). Thus given a MNS 
M represented as M = {Si, S2, ■ ■ ■ , S r } we will first take the residues with 
respect to the roots 0(lk), for sets Sk of cardinality 2, in arbitrary order, 
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[1,2,3,4,5,6,7] 



[10,11] 
[9, 10, 11] 
[8,9,10,11] 



[1,2,3,4,5,6,7,8,9,10,11] 
Figure 8. Irreducible components of a MPNS in A\q 




Figure 9. Order of nodes in the tree represented in Figure 8, 
according to Proposition 6.5 



then with respect to irreducible sets associated to sets of cardinality 3, 
etc. The procedure of ordering roots coming from a MNS M = {Sj} ac- 
cording to the cardinality of the set Sk is called OrderThetas. Furthermore 
we will at the same time keep track of the order of the pole for calcu- 
lating an iterated residue of a function (f> = P/ YiiLi a i m the- procedure 
FormalPathAwithOrders. 



7. Volume and partition function for the system A n _i 

7.1. The formulae to be implemented. Let E be an n-dimensional vec- 
tor space with basis ei (i = 1, . . . , n) and consider the set 

K-n = {&i - ej | 1 < i < j < n}. 

These are the positive roots for a system of type A n -\. The number of 
elements in K, n is N = n(n — l)/2. Note that K, n is also the set of vectors 
in a complete graph with n nodes. 
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We let V be the vector space generated by the elements in )C n . Then V 
has dimension n — 1 and it is defined by: 



n 



i=l 



In our procedures, a vector v of length n such that Yli=i v i = ® wm be called 
an A-vector and written as v = [v±,V2, ■ ■ ■ ,v n ]. The lattice spanned by fC n 
is simply 



n n 



V z =lh = Y J h t e t € 7L n | ^ = L 
I i=i i=i J 

It is well known and easy to prove that /C n is unimodular. The cone C(/C n ) 
generated by /C n is simplicial with generators the n — 1 simple roots ei — 
62, e 2 — e3, . . . , e n -i — e n . This cone is described as: 

C(/C n ) = {^4- vector u = [ui, «2, . . . , v n ] \ vi + V2 + ■ ■ ■ + Vi > for all i}. 

Keep in mind that our vector v satisfies the condition 

vi+v 2 -\ h v n -i + v n = 0. 

We choose on V the measure dh determined by V%. Let v be in the cone 
C(/C n ). We are interested to compute the volume volz,K n ( v ) °f the polytope 

If h is a point in V with integral coordinates then we are also interested in 
computing the number N^ n {h) of integral points in ILjc n (h). 

We apply the formulae of Theorem 3.3. Since KL n is unimodular, the set 
F can be taken as F := {0} (Remark 3.8). 

Since V is contained in E, then we have a canonical map E* — ► V* given 
by restriction. Define U = V* as in the general setting. We identify U with 
IR n_1 by sending u £ W 1 ^ 1 to u = J2i=i u i e% e E* ■> where e* is the dual basis 
to ej. Thus the root — ej (1 < « < j < n) produces the linear function 
Ui — Uj on U, while the root e$ — e n produces the linear function U{. 

Definition 7.1. Let u = X^I=i v « e « S V be a vector with real coordinates. 
Let h = Y17=i hi e i ^ V be a vector with integral coordinates. Then for 
u <E U define: 



Ja(w)(u) 



FA(h)( U) - _ x 



rnu 1 u * rii<i<j<n-i(^ - «j) 
nr=i 1 (i+^) /ii+n_1 ^ 



nr=i ^rii<i<j<n-i(^-^) 

Theorem 7.2. Lei t be a chamber ofC(K, n ). 



30 



M.W. BALDONI, M. BECK, C. COCHET, M. VERGNE 



• For v G c, we have 

^oh,ic n (v) = JK c (J A (v)). 

• For ft 6 Z"nc, we have 

N Kn (h) = JK c (F A (h)). 

Proof. The first assertion is the general formula. 

The function F(0, h){u) = e^ / H aeA (l ~ e" (a ' M> ) for the system K n is 

e Y,i=l u i v i 



F(0, h)(u) = 1 

iETi'a-e-^jni^n-ia-c-^-^) 

Note that the change of variable 1 + Zi = e Ui preserves the hyperplanes 
Ui = and Ui = Uj. After the change of variable, we get 

(11) F(0, h)(u) = 1W V + _ wn »-i ' 

But Zi = e Ui — 1 leads to dzi = e Ui du% = (1 + Zi)dui and hence we obtain the 
desired exponent hi + n — i — 1 thanks to the formula involving Jacobians 
in Proposition 3.2. □ 

In order to implement these formulae, we first have to describe the set 
V(v,)C n ) (Section 7.2), then calculate the iterated residue formulae associ- 
ated to these paths (Section 7.3). Below we explain how these computations 
fit together to get a global procedure for the Kostant partition function for 
A n _i (Section 7.4). As a short digression, we will explain how we adapted 
our program to deal with formal parameters (Section 7.5). 

7.2. The search for maximal proper nested sets adapted to a vec- 
tor. We now look for maximal proper nested sets adapted to a vector fol- 
lowing the general method as outlined in Figure 7: we will begin by listing 
all possible /C n -admissible hyperplanes. The usual height function is 

n-l 

ht(v) = y](n — i)vi 
i=i 

which takes the value 1 on all the simple roots, and hence the value j — i on 
ei — ej. We deform ht slightly in order to have a function taking different 
values on all roots: If two elements ej — ej and — eg are such that j — i = 
f - ft, we decide that ht(ej — ej) < ht(e^ — eg) if i < k. 

If P is a proper subset of {1,2,..., n} and v is an ^4-vector, we denote by 
(up,v) the linear form XlieP^' an( ^ by Hp the hyperplane 

Hp := {v e V\(u P ,v) = 0}. 

We will see shortly that all /C n -admissible hyperplanes are obtained in 
this way, that is giving a proper subset P of {1, . . . ,n}. Observe that the 
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hyperplane Hp is equal to the hyperplane Hq determined by the complement 
Q of P. We denote 

K{P) := {e; - ej | 1 < i < j < n ; i, j G P} C K n . 

Note that JC(P) is the positive system A\ P \_ 1 , where the positivity is induced 
by the lexicographic order. 

Lemma 7.3. • The hyperplane Hp is a JC n - admissible hyperplane. 

• The set K n n Hp is the union of JC(P) and )C(Q), where Q is the 
complement of P in {1, 2, . . . , n}. 

• Every JC n - admissible hyperplane is of this form. 

Proof. The first two assumptions are easy to see. We prove the third by 
induction on n, the case n = 2 being trivial. Let H be a /C„-admissible 
hyperplane. Let a be a root in H. Renumbering the roots, we may assume 
that a = e n - \ — e n . The map q sending ei to if i < n and e n to e„_i sends 
the set K n \ {a} to /C n _i. The space H/Ra becomes a /C„_i-admissible 
hyperplane. It is thus determined by a subset P' of {1,2, ... , n — 1}. If 
P' does not contain n — 1, the hyperplane H is equal to the hyperplane 
determined by the subset P' of {1, 2, . . . , n— 1}. If P' does contain n— 1, then 
the hyperplane is equal to the hyperplane determined by P = P'u{n}. □ 

We now proceed to the detailed description of our algorithm. Recall our 
description of an ^4- vector as an array v = [vi,V2, ■ ■ ■ ,v n ] with J27=i v i = ®- 
Referring to Figure 7 we need to check if the vector is in the cone C(/C n ), 
that is, Y^rj=i v j — f° r 1 < i < n — 1. This is done by using the procedure 
CheckVector (v) , which gives an answer true or false. 

For the system JC n the highest root 9 is equal to 

9 = [1,0,... ,0,-1] G R n 

and computed with the procedure theta. 

At this point, we need to list all hyperplanes that are in Hyp(f,/C n ). 
This is done in the procedure TwoSets(v), that we are about to describe. 
As explained in Lemma 7.3, each hyperplane is determined by an equation 
Yliei a ' 1 = 0- I* therefore produces a set of two lists P, Q, where P = [i € I] 
and Q = [i $ I \ 1 < i < n]. Note that P and Q are sorted. To verify that 
such a hyperplane is in Hyp(f,/C„), we need to test if (up, 8) is not zero (0 
is not in the hyperplane) and if (up,v) x (up, 6) is non-negative (0 and v 
are on the same side of the hyperplane). 

Furthermore, the procedure ProjH(v, H) constructs the vector 

• i \ ( U P, V ) a 

proj^(w) = - 



(up, e) 

that we represent as {[ui,P], [v2,Q]}- Each of the vectors v±, V2 is an A- 
vector (sum of coordinates equal to zero). So the last condition for H being 
in Hyp(w,/C n ) is that v 1 £ C(K(P)) and v 2 € C(/C(Q)). 
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Hence a hyperplane H is in Hyp(v,/C n ) if it satisfies the series of condi- 
tions: 

(up, 9) / with Hvalue(theta(n),P) / 0, 
(up,v) x (up, 9} > with CheckSide(v,P) = true, 
v± G C(/C|p|_i) with CheckVector(vi) = true, 
V2 G C(/C|Q|_!) with CheckVector(v 2 ) = true. 



The procedure CheckList(v, H) implements all these sub-routines. It is used 
in the procedure TwoSets(v), computing all elements of Hyp(w,/C n ). We 
combine TwoSets with a procedure named TwoVector to finally get the 
procedure TwoVectors(v) determining all hyperplanes in Hyp(v,/C n ) and 
projections of v on these hyperplanes. 

We now have to perform the next step of our algorithm. Let { [v i , K\] , [v2 , K2] } 
be the output of TwoVectors(v). Then we construct the MNSs for 
and [v2, K2], and go on recursively until the procedure stops. These iterated 
steps are done by the procedure Splits. 

Finally the procedure MNSs(v), computing all MNSs for a given vector 
v, works as follows. We begin by building the first seed of MNSs with the 
procedure MNS1, containing the regularization of the result of TwoVectors. 
We then call repeatedly the procedure AllNewMNSs, which performs the 
regularization of the output of Splits. 

7.3. Residues associated to maximal proper nested sets. An element 

M in V(v,K n ) is represented as a collection M = {Ki, K2, . . . , K n -i} of 

(n — 1) subsets of [1,2,... ,n\. As we have said in Remark 6.6, given a 

maximal proper nested set M := {K\, K2, ■ ■ ■ , K n -i} we associate to it an 
> 

ordered basis 9(M) of V (procedure OrderThetas). If p = [cti, ct2, ■ ■ ■ , at n -i] 
is the list of roots singled out by our procedure, then a.\ is an element 
associated to a set Ki of cardinality 2 and a n _i = 9. We identify the root 

— e n to the linear function z\ on C™" 1 and the root — ej to z\ — Zj. 

Let h be a yl-vector with integral coordinates. Let us consider the Kostant 
function 

lll<i<i<n-lV z * z 3> X lli=l z * 

(Definition 7.1). To compute Njc n (h), we will have to compute 
resM</> := resQ, n _ 1= ores a2= o • • • res ai= o4> 

with 4> = J^A{h). Using Proposition 6.5, we know in advance the order of the 
pole in «fc = of the function obtained after taking the first [k — 1) residues. 
These orders are recorded in the procedure FormalPathAwithOrders. 

If ot\ = Zi — Zj, we can replace — after taking the residue at Zj = Zj — the 
variable Zj by the variable Zj in all the other roots. Thus we get rid of the 



PARTITION FUNCTIONS FOR CLASSICAL ROOT SYSTEMS 33 

variable z%. The procedure NewR produces the ordered path resulting from 
all these substitutions. 

Recursively, we will have to compute the residue at Zi — Zj = of an 
expression 

(12) / = A{Zi,i£ ^\ m , 

Yli,jeL;i<j( Z i ~ Zj) mi,J T\i£L Z i 

where L is a list of indices taken in {1, ... ,n — 1}. Denote by maxi the order 
nii (h j of the root Zi — Zj (the exponent maxi is recorded in the procedure 
FormalPathAwithOrders). Note that computing the residue is exactly the 
same as computing the coefficient of z of degree maxi — 1 of the expansion 
of / x (zi — Zj ) maxt at Zi = z + Zj . Let us describe in detail the procedure 
ComputeRes, performing this task. 

For j G L\ {i ,io}, let et ej = m i(hj if i < j and ej = m jtio if i > j. 
Then / x (z io — Zj ) maxt can be written as g x h x R where 

(13) g = = A ^^ L l x 1 _ 

( 14 ) h = Y\ 7 \T x 4o~> 

(15) R = (-rpi>i e i. 

Let Zi = z + Zj . So to get the desired residue, we need to calculate 
the expansion of g and h at z = 0. More precisely if g = Yli^Q 1-1 9i z% an d 
h = ^^o" - hiZ 1 , then the coefficient of degree maxi — 1 of fx (zi — Zj ) maxt 
is simply R x ^I^T" 1 di^maxi-i-i- Let us describe how the procedure 
ComputeRes performs this task. 

Rewrite the fraction h defined in Equation (14) as B x h, where 



B 



njeL;j^i ,j ( Z 3 Z jo) ej z jo 
~ h = 1 1 

Consequently, to expand h as a function of z, we only need to expand h. 
This is done in the procedure Coef f Bin using the binomial coefficients. In 
the procedure Coef f Fun we calculate the expansion at Zi — Zj = of the 
fraction 



(16) fx(z H) -z Jo r™ xzg>x (*o-*jT j 

= g x Rx (-1)^*0 jo 
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Finally the procedure ComputeRes performs the sum over i ranging from 
to maxi of 

(S ■= (-l)^o^) x B 
x (the component of degree i of Coef f Fun) 
x (the component of degree maxi — 1 — i of Coef f Bin). 
Rewrite this as the sum over i of 
R 

x (the component of degree i of Coef f Fun) x R x S 

x B x (the component of degree maxi — 1 — i of Coef f Bin), 

or, equivalently, as ^ImT -1 x 9i x h max i-i-f- this is exactly the desired 
coefficient. 

Remark 7.4. For residues along roots of type Zi instead of Zi — Zj the 
procedure ComputeRes also calls procedures srCoeffFun and srCoeffBin, 
similar to CoeffFun and CoeffBin. 

7.4. The procedure MNS_KostantA. We finish the section dedicated to 
A n -\ by giving the global outline of the procedure MNS_KostantA(v) com- 
puting the Kostant partition number of a vector v lying in the root lattice. 
We begin by slightly deforming v so that it lies on no admissible hyper- 
planes, with the command v' := Def Vector(v, n). We compute all MPNSs 
for v' with the procedure MNSs(v'). 

Given such a MPNS M = {Sk}, we extract the highest roots of its ir- 
reducible components with the call R := ThetaMNS(M). We obtain a set R 
where each element of R is a root represented as [i, j] together with the cardi- 
nality of the set Sk it comes from. We then transform this set R into a path p 
keeping track of order of poles by setting p := FormalPathAwithOrders(R). 

Finally we compute the residue associated to this path with the command 
OneIteratedResidue(p, v,n). Summing all these residues over the set of 
MNSs, we obtain thanks to Theorem 7.2, the desired partition number for 
v. 

Let us describe in detail the procedure OneIteratedResidue(p, v, n) com- 
puting the iterated residue along a path p for a vector v lying in the root 
lattice for A n -\. We first compute the Kostant fraction (second item of Def- 
inition 7.1, procedure KostantFunctionA). Then we replace in the path all 
roots Zi — z n by z,- L (with Kpath). We also build a upper bound for the orders 
of the roots rriij (with Multiplicity). Keep in mind that the exponent 
maxi needed in the residue calculation is computed a priori, and our com- 
putation seems quite optimal. Then we compute iteratively the residues, 
using the procedure ComputeRes (Section 7.3). Note that at each step we 
have to update the list of orders (with MultRoots) and the list of remaining 
variables (with ListOf Variables). 

7.5. Parametrized version of the algorithm. Our algorithm can work 
with formal parameters, only needing slight modifications of procedures. We 
are then able to compute directly the polynomial h i-> N/c n (h) giving the 
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number of integral points in the polytope Hjc n (h), on the chamber deter- 
mined by h (this chamber is easily computed). As a consequence we can 
easily get the Ehrhart polynomial t i-> Njc„(thi, . . . ,th n ) of the polytope. 
See [13]. 

Now let us outline how this modified program works. Given an ele- 
ment h = (h±, . . . ,h n ) of the root lattice for A n -±, we want to compute 
the Kostant partition function for the vector (hi, . . . ,h n ), when h varies in 

z n nc . 

Recall that we have to perform the residue at Zi — Zj = of the fraction 
defined in Equations (11) and (12). Note that the numerator of the fraction 
J~A(h) contains terms of the form (1 + Zj ^ h i+ n - 1 - 1 that we must formally 
expand. For any parameter b, we write at Zi = z + Zj , 

(maxi-l , \i\ 

1+ E (j) {tTz-J j +°(^)- 

Hence, for any parameter b, the procedure CoeffFun (see Equation (16)) 
now computes the expansion at Zi — Zj = of the fraction 



/ >< k°x n (*»- 



maxi—l /, 





maxi — 1 

JO J 



For the residue along a root of type Zi instead of Zi — Zj , the procedure 
srCoef fFun has been modified in a similar way. 



8. The type B n 

8.1. The formulae to be implemented. Consider a vector space V with 
basis ei, e2, • • • , e n . We choose on F the standard Lebesgue measure dh. 
Let 



B n = {e i \l<i<n}U{e i -e j \l<i<j<n}U{e i + e j \l<i<j< n}. 

Then B n is a positive roots system of type B n and generates V. The number 
of elements in B n is iV = n 2 . We denote by C/ the dual of V. The lattice V% 
generated by roots is equal to Z n , so the constant vol (V/Vz, d/i) = 1. 

The cone C(B n ) is simplicial and spanned by the n simple roots e\ — ei, 
e 2 - e 3 , . . . , e n _i - e n , e n . A vector v = [vi, v 2 , ■ ■ ■ , v n ] is in C(H n ) if and 
only if it satisfies the inequations v\ + ■ ■ ■ + Vi > for alH = 1, . . . , n. 

Let v be in the cone C(B n ). Consider the polytope 
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n B„(«) = < (x a )a > ^ x a a = V 



If h is a point in V with integral coordinates, we are interested in computing 
the number Ng n (h) of integral points in 



Let Uz be the lattice dual to Vz- We identify the torus T = U/Uz = 
W l /Z n to (S* 1 )™ by 



If G is a representative of g = (gi, g2, ■ ■ ■ , <? n ) G 7\ and /i = Y17=i ^* e * ^ n ^> 
then e (h,2nV^iG) j g e q Ua i ^ Q rj™ =1 g^ = As the set B n is not unimodular, 
the sets T(a) are not reduced to 1. 

Example 8.1. Let a be the basic set {e\ + e2,e\ — e{\ for F>2. Then 



We now determine a set F containing all sets T(a). 

Lemma 8.2. Let a be a basic subset of B n . Assume g G T(a). Then all 
the coordinates of g are equal to ±1. Furthermore, if g is not 1, there are 
at least two coordinates of g which are equal to —1. 

Proof. We prove this by induction on n. For £>2, we have seen this by direct 
computation. 

Let a be a basic subset of B n . Assume first that a contains a root e*. 
Up to renumbering, we may assume that this root is e n . Then the basis a 
produces a basis a' of B n -\ by putting e n = 0. Let g = (gi,g2, ■ ■ ■ ,g n ) i n 
T(a). We see that g' = (gi, g2, ■ ■ ■ , g n -i) is i n T(a'). Thus, by induction 
the first n — 1 coordinates of g' are equal to ±1. But since e n is in a we get 
1 = g n . Note that g ^ 1 if and only if g' ^ 1, hence by induction hypothesis 

has at least two coordinates not equal to 1. 

Consider now the case where a does not contain any root ej. Up to 
renumbering, it contains a root e„_i — e n or e„_i + e n . 

Let us examine first the case where a contains the root a = e n _i — e n . 
Let g = (gi,g2, ■ ■ ■ ,g n -i,9n) in T(a). This implies g n ^i = g n . Consider the 
map q sending to if i < n and e n to e n -±. Then q sends a \ {e„_i — e n } 
to a basis <r' of B n -\. The element </' = (gi,g2, ■ ■ ■ ,9n-i) is easily seen to 
belong to T(a'). Indeed if a equals ej ± ej with 1 < i < j < n, this is by 
definition. On the other hand q{ei ± e n ) = e^ ± e„_i and g n -i = imply 
that gig^-i coincides with the value of gig^ 1 ■ By induction hypothesis, all 
coordinates of g' are equal to ±1. Moreover g ^ 1 if and only if g' ^ 1, so 
that <7 is of the desired form. 

Finally, the same argument works if a contains a = e n _i + e n , by consid- 
ering the map q sending ej to e^ if i < n, and e n to — e n _i. □ 




2» = {(1,1), (-1,-1)}. 
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Definition 8.3. If / is a subset of {1,2,..., n} with at least two elements, 
we consider the set F(I) := {(51,52, ■ ■ ■ ,9n) \ 9i = -M G I;gj = l,j £ I}. 

We define F C T to be the finite subset of T union of such sets F(I) 
together with the identity (1, 1, . . . , 1). 

Let v = Y17=i Viei G V be a vector with real coordinates and h = 
Yli=i hiei € V a vector with integral coordinates. We will compute the 
normalized volume of IIg n (v) and the number of integral points in ITg n (h) 
using Theorem 3.3. Thus we introduce the function Jb(v) on U defined by: 

En 
i= i 

Jb(v)(u) 



n?=i u * rii<i<i<n(^ - u i) rii<i<j<n(^ + n i) ' 

For 5 = (51, 52, • • • ,<?n) G F and h E VzD C(B n ) the Kostant fraction (3) is 
the function on [/ defined by: 



F B (g,h)(u) = 



We have then 



1 



Theorem 8.4. Lei c be a chamber ofC(B n ). 

• For any v £c, we have 

yoh,B n {v) = JK C {Jb{v)) ■ 

• For any h £ Vzflc, i/ie uaZtie 0/ i/te partition function is given by: 

N Bn (h) = J2^MFB(9,h)). 

g&F 

As in the case of A n , we will use the change of variable 1 + Z{ = e Ui to 
compute more easily Ns n {h). However, let us note that this transformation 
does not leave the hyperplane Uj + uj = fixed. This hypersurface is trans- 
formed into the hypersurface Zi + Zj + ZiZj = 0. So we use the expression of 
JK C as an integral over the cycle H(c) defined in Theorem 4.14. This cycle 
(its homology class) is stable by the transformation e Ui = 1 + Zi which is 
close to the identity. Thus define the following function on U: 



FB(g,h)(z) = 



IE = i(i + ^) hi+aB - i - 1 xIE=i^ 



EKUt 1 + Z i ~ 9i) X ]ll<i<j<n( 1 + z i ~ 9i9j{ 1 + z j)) 
1 



]ll<i<i<n( 1 + z i)( l + z i) - gigj ' 

Performing the change of variables e Ui = 1 + Z{ on the function Fs(g, h){u) 
and computing the Jacobian, Theorem 3.3 becomes: 
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Theorem 8.5. Let c be a chamber ofC(B n ). 

• For any v G c, we have 

voh,B n (v) = JK C {Jb{v)) ■ 

• For any h^VzHc, the value of the partition function is given by: 

NB n (h) = Et^TvT / FB(g,h)(z)dz. 
^ F (2vrV-l) n Jh(c) 

As in the case of type A, in order to implement these formulae we first 
have to describe the set V(v,B n ) (Section 8.2), then we will explain how 
the integral over the cycle H (c) is calculated similarly to an iterated residue 
formula associated to these paths (Section 8.3), using an estimate of the 
order of poles. Finally we explain how these computations fit together to 
get a global procedure for Kostant partition function for B n (Section 8.4). 

8.2. The search for maximal proper nested sets. A height function is 

n 

ht(v) = ^(n + 1 - i)vi 
i=i 

which takes value 1 on all simple roots. We will deform ht later on in order 
to have a function taking different values on roots. 

We now proceed to describe hyperplanes for B n . If P = [P + ,P~] are 
two disjoints subsets of {1, 2, . . . , n}, we denote by (up, v) the linear form 
X^eP+ Vi ~ YljeP- v r Consider the hyperplane 

Hp ={«€ V, (u P ,v) = 0} 

in V. It is equal to the hyperplane determined by the reverse list [P~,P + ]. 
Thus to each set P = {P + ,P~} of two disjoint sets P + , P~ such that at 
least one is non empty, we associate a hyperplane Hp. 

We denote by Z the complement of P + U P~ in {1, 2, . . . , n} and by B(Z) 
the subset of B n defined by 

B(Z) = { ei \ie Z}U{ei±ej\l<i<j <n;i,j G Z}. 

This is the positive root system B>i z i, with the positivity induced by the 
lexicographic order. 

Let 1C(P + ,P~) be the subset of B n defined by 

{&i - ej 1 1 < i < j < n ; i, j € P + } 
U {ei + e k | i G P + , k G P'} 
U {e k - e e | 1 < k < £ < n; k,£ G P~}. 

Note that by defining /j = if i G P + and fk = — e|p-|-fc+i if fe G P~ , the 
set )C(P + ,P~) coincides with 
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{fi- fj\l<i < j <n; i,j£P + } 

u {fi-f k \ieP + ,keP-} 

U {fk-fe\l<k<£<n;k,£eP-}. 

Thus the set )C(P + ,P~) is a positive root system of type ^4|p+| + |p-|_!. 
However the positivity is induced by the lexicographic order on P + and the 
reverse lexicographic order on P~ . Observe also that Hp is the vector space 
spanned by /C(P+, P~) U B(Z). 

Lemma 8.6. • The hyperplane Hp is a B n - admissible hyperplane. 

• The set B n n Hp is the union of B(Z) and K.(P + ,P~). 

• Every B n - admissible hyperplane is of this form. 

Proof. The first two assumptions are easy to see. We prove the third as- 
sumption by induction on n, the case n = 2 being trivial. Let H be a 
^-admissible hyperplane and let a be a root in H. There are 3 possibil- 
ities for a: up to renumbering roots, we can consider the cases a = e n , 
a = e n -i - e n and a = e„_i + e n . 

In the first case, the map q sending to ej if i < n and e n to maps 
the set B n \ {a} to B n -\. The space H/M>a becomes a ,6 n _i-admissible 
hyperplane. It is thus determined by P' = [P + , P ~], where P + and P ~ 
are two disjoint sets contained in{l,2,...,n — 1}. Then the hyperplane H 
is equal to the hyperplane determined by [P' + , P'~]. 

In the second case, the map q sending e« to ej if i < n and e n to e n _i 
sends the set B n \{a} to B n -\. The space H/Ra becomes a £> n _i-admissible 
hyperplane. It is thus determined by P' = [P + , P ~\. If neither P + nor 
P _ contain n — 1, the hyperplane H is equal to the hyperplane determined 
by [P + , P ~\. Otherwise assume that for example P + contains n—1. Then 
the hyperplane H is equal to the hyperplane determined by [P + , P~], where 
P+ = p'+ u {n} and P~ = P'~. 

In the third case, the map q sending ei to 6^ if % < rx and e n to — e n —\ 
sends the set B n \{a} to B n -\. The space HfRa becomes a S Tt _i-admissible 
hyperplane. It is thus determined by P' = [P + , P ~\. If neither P + nor 
P ~ contains n — the hyperplane H is equal to the hyperplane determined 
by [P + , P ~]. Assume that P + contains n — 1. Then the hyperplane H is 
equal to the the hyperplane determined by [P + ,P~], where P + = P + and 
p- = p'~ u {n}. □ 

We now give a detailed description of our algorithm computing maximal 
nested sets. We describe a vector as an array v = [vi, V2, ■ ■ ■ , v n ]. To check 
if v is in the cone C(B n ), we need to verify if J2)=i v j > for 1 < i < n. 
This is done by the procedure CheckBvector, which returns the answer true 
or false. 

For the system B n the highest root 9 B {n) is equal to 
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B (n) = [l,l,O,O,O,...,O]. 

We recall here that P is divided in two sets P + U P~ , one of them being 
non empty. The first task is to list the hyperplanes in ~H.yp(v , B n ) . This set 
of hyperplanes is obtained by the command line AllPossibleBwalls(v). 
The input of this procedure is the vector v. The output is a set of elements 
P = {P + , P~}, where P + = [h,i 2 , ■■■i p ] and P~ = [ji,j 2 , ■ ■ ■ ,j g ] are two 
ordered disjoint lists made from indices taken in {1, . . . ,n}, with at least 
one of P + or P~ being non empty. Let (up, v) = YlieP+ v i ~ Sjgp- v j be 
the normal vector to Hp. Then as stated in Lemma 5.2 we need to test if 
(up,9 B (n)} is not zero and if (up,v) x (up,6 B (n)} is non negative. 

We then construct the vector 

■ / A (UP,V) B 

pmiM=v - ( Up ,eB(n)) e {n) - 

This vector is represented as {[v\,P + ], [v2,P~], [w,Z]}. The sum of coor- 
dinates of vi is equal to the sum of the coordinates of v 2 . Now Z is the 
ordered list \k\, k 2 , ■ ■ ■ , kg] of complementary indices to P + U P~ and 

w = \pToj H (v)[ki], . . . ,pioj H (v)[ke]]. 

Note that the equations of the cone C()C(P + , P~)) can be given in the 
convenient form v\ © v 2 G C(IC(P + , P~)) if and only if CheckBvector(vi) 
and CheckBvector(v 2 ) are true. Equations of the cone C(B(Z)) are given 
in the form w G C(B(Z)) if and only if CheckBvector(w) is true. 

Thus the condition that H is in Hyp(i;,£> n ) is equivalent to the series of 
conditions: 



(u P ,e B (n)) + 0, 

(u P ,v) x (u P ,e B (n)} > 0, 

CheckBvector(vi) = true, 

CheckBvector(v 2 ) = true, 

CheckBvector(w) = true. 

Those five conditions are checked by the command line CheckBwall(v, H), 
that gives an answer true or false. 

Remark 8.7. We can first construct all disjoint subsets P + , P~ of {1,2,..., n} 
and test these five conditions successively on all of them. However it is 
highly desirable to throw away a priori a great number of these partitions 
by noticing the following restrictive conditions on the possible lists to be 
considered. 

Let {P + ,P-} = {[i i,i 2 , . . . i p ], [ji,j 2 , ■ ■ ■ ,jq]} be a set of two disjoint 
subset of {1,2, ...n} represented as lists with strictly increasing indices. 
Let Z = [hi, k 2 . . . , kg] be the list of complementary indices to P + U P~ in 
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{1, . . . , n}. The following linear forms are positive on the cone C(JC(P + , P )) 
generated by )C(P + ) and K.(P~): 

+ vi s > for all 1 < s < p, 
+ v jt > for all 1 < t < q, 
+ v ks > for all 1 < s < I. 

Note that all the above linear forms take positive values on 6 B (n). We 
employ Lemma 8.6. Thus if v[ii] < 0, the index i\ cannot start the list P + 
of an element {P + ,P~} in AllPossibleBwalls(v) and we reject all such 
{P + ,P-}. 

Similarly assume that we have constructed a list of indices 12] satisfying 
conditions v[i\\ > and v[i\\ + v[i2\ > 0. Then if v[i\] + v[i2\ + v[is] < 0, a 
list starting with 12, 13] cannot be the first three indices of the component 
P + of an element {P + ,P~} in the set AllPossibleBwalls(v) and we skip 
it right away. 

This achieves the description of the procedure AllPossibleBwalls. We 
now have to perform the next step of our algorithm. As for type A we build 
MNSs iteratively. At each step we get a set of partial MNSs, to which we will 
apply recursively our algorithm. Note that after Lemma 8.6 the intersection 
of a ^-admissible hyperplane Hp with B n is the union of a system of type 
A and a system of type B. 

The part of the MNS coming from the subsystem of type A is computed 
with the procedure AddAnests. It performs a reordering of the result of a 
call to the procedure MNSs described in Section 7.2. 

The part of the MNS coming from the subsystem of type B is com- 
puted with the procedure Bsplits, calling the previously described proce- 
dure AllPossibleBwalls. 

Procedures AddAnests and Bsplits are enclosed in MoreNSs, thus giving 
a new iteration of the process. After regularization of the result we hence 
get a procedure named AllNewNSs, performing a new step in the building 
of MNSs. 

Finally the procedure B_MNSs, computing MNSs for a given vector v for 
type B, is the following. First, we use a procedure named B_NS1 to calculate 
the first seed of all MNSs. After, repeated calls to the procedure AllNewNSs 
build the desired MNSs. 

8.3. Residues associated to maximal proper nested sets. A proper 
maximal nested set M gives rise to an ordered basis Oj, and a cycle H(M). 
We need to compute 



Vn + Vi 2 H 

v h + v h + ■ ■ ■ 
Vki +v k o-\ 



I F B {g,h)(z)dz 

H(M) 
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where 

H(M) := {z,\(a l ,z)\=e l }. 
The function z Ts^g-, h)(z) is deduced from the function Fs(g, h)(u) in 
the space -R4 by the change of variable e Ui = 1+Z{. Thus its denominator is a 
product of factors, either of the form Z{ corresponding to the root Ui, or of the 
form Zi — Zj corresponding to the root Ui — Uj or Zi + Zj+ZiZj corresponding to 
the root Ui + Uj. We denote by u(z) the point with coordinates U{ satisfying 
e u i = 1 + Zi . 

We start integrating our function Tb{9i h)(z) over the smaller circle | (a\, z) \ = 
ei keeping the other variables fixed. By our condition on the cycle, the 
function we integrate has poles on the domain |(ai,z)| < e\ only when 
a\{u(z)) = 0. If oti{u(z)) = ui — Uj or a\{u(z)) = Ui, the poles are ob- 
tained for Zi = Zj or zi = 0. If a\(u(z)) = ui + Uj, the pole on the domain 
|(ai,z)| < ei is obtained for Zi = —Zj/(1 + Zj). Thus we compute the in- 
tegral over the circle by the residue theorem in one variable, and proceed. 
From the general theory, the poles of the function we obtain, replacing Zi 
by one of the values above are again of the same form with respect to the 
remaining variables, as is easily checked. 

As in case A n —i, for a root a = Ui (resp. a = Ui =t v,j) we can replace 
after taking the residue at a = the variable Zi by (resp. by TZj) in all 
other roots. Thus we get rid of the variable Zj. The procedure FormalPathB 
produces the ordered path resulting from all these substitutions. 

In the case of type B we compute the residue by directly checking the order 
of the pole at a = 0, and then using differentiation. The program works in 
the same way with parameters. The function obtained is locally polynomial 
with polynomial coefficients depending of the parity of the integers hi. 

8.4. The procedure MNS_KostantB. We finish the section dedicated to 
B n by giving the global outline of the procedure MNS_KostantB(v) comput- 
ing the Kostant partition number of a vector v lying in the root lattice of 
B n . We begin by slightly deforming v so that it lies on no wall, by set- 
ting v' := Def VectorB(v, n). We then compute all MNSs for v' with the call 
B_MNSs(v') (Section 8.2). For every MNS M, we extract the list R of highest 
roots of its irreducible components by setting R := BthetaMNS(M). We sort 
these roots by their height with the command line R' := BorderThetas(R, n). 
We then transform the list of roots R' into a path p by setting p := FormalPathB 

Now remark that our procedures are designed to take residues along pos- 
itive roots, using the fact that res_ a = — res a for any root a. The sign that 
appears (more precisely —1 to the power the number of negative roots in 
the path p) is computed with the procedure PathSign(p, n). 

Then for every g in F we do the following. The iterated residue along the 
pathp and for g is obtained by the command line OneIteratedBresidue(p, g, v, 
Let us briefly describe its implementation. We first compute the Kostant 
fraction (second item of Definition 7.1, procedure KostantFunctionB). Then 
for every root of the path we apply the procedure ComputeOneResidue 
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(Section 8.3) and update the order of the pole with a procedure named 
OrderPoleB. 

Finally summing all products PathSign(p, n)xOneIteratedBresidue(p, g, v, n) 

over the sets of g's and of M's, we get the desired result. 

Remark 8.8. Let us fix a list R' = [a>i, . . . ,a n ] of ordered roots coming 
from a MNS, and an element g. We say that R' and g are compatible if the 
following condition is satisfied. If indices of monomial(s) of have not yet 
occured among indices of roots ag with £ < k, then g must satisfy g ak = 1 
(that is g%gf l = 1 if afc = ej ± ej and gi = 1 if = e^). Note that the 
iterated residue for g and for the path p associated to R' is zero if g and 
R' are not compatible. Hence summing only over g's that are compatible 
with a given list R' saves useless computations. The check of compatibility 
is performed by the procedure ListAndGAreCompatible(R', g, n). 



Consider a vector space V with basis e±, e2, . . . , e n . We choose on V the 
standard Lebesgue measure dh. Let 

C n = {2ei 1 1 < i < n} U {ej — ej 1 1 < i < j < n} U {ej + ej 1 1 < i < j < n}. 

Then C n is a positive roots system of type C n , and generates V. The number 
of elements in C n is N = n 2 . Note that elements of C n and B n are propor- 
tional, so they determine the same hyperplane arrangement and the same 
chambers. 

Let L be the lattice defined by Zei © Ze2 © • • • © Ze n . We remark that 
the lattice generated by C n is the sublattice of index 2 in L consisting of 
all elements v = [v± , v-i , . . . , v n ] with integral coordinates and such that the 
sum Y^l=i Vi 1S an even integer. A Z-basis of Vz is, for example, 

Z(ei - e n ) © Z(e 2 - e n ) © • • • © Z(e„_i - e n ) © Z(2e„) , 
so vol ty/Vz) = 2. 

The dual lattice U% is the lattice of vectors 7 = (71,72, ■■■ ,7n) such 
that 7j are half integers and such that 7, + 7^ is an integer for all i, j. 
The set C/^/Zei © • • • © Ze„ is of cardinality 2 with representative elements 



(0,0,... ,0,0) and (1/2, . . . , 1/2). 

As before, we identify the torus f = U/(Zei © • • • © Ze n ) = R n /Z n with 
(5 1 )™ by 



T = T/{±1} = [//%. 

Let G be a representative of g = (51,52, • • • , 9n) £ T and /i = Yli=i ^« e « 
in F z . Then e (h^^lG) ig equal to J]™ =1 ^ ! = 5 \ This function is well 

defined on T = T/{±1} since X^r=i ^» ^ s even - 
For a a basic subset of C n , define 



9. The type C, 



n 




Then 
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f{a) = [g^f e ( a ^^iG) = 1 for allQ , G a j _ 

As the set C n is not unimodular, sets f(<r) are not reduced to 1. 

Lemma 9.1. Let a be a basic subset of C n . Then T(a) C {±l} n . 

Proof. We prove by induction on n that if a is basic then the condition 
g = (<7i, . . . ,g n ) G T(c) forces g| = 1 (1 < i < n). In other words g a = 1 
for all long roots a. If so then gi = ±1 for all i. The base of the induction, 
that is C2, is straightforward and we omit it. We thus proceed considering 
various possibilities for our a. 

If there exists a long root in a we may assume that this long root is 2e n . 
We embed the system C n _i in C n via the first (n — 1) coordinates. Then 
the basis cr of C n produces a basis cr' of C n _i consisting of roots {ej ± ej G 
cr 1 1 < i < j < n — 1}, of roots {2ej £a\l<i<n — 1}, and of roots 
{2ej I ej±e n G a;i / n}. It is easy to see that the elements (g±, g2, ■ ■ ■ , g n -i) 
are in T(a J ). Indeed gf = 1 if i e„ G <r as gigt 1 = 1 an d 5n = lj an d 
similarly g| = 1 if 2ej G cr. Thus by induction we obtain gf = 1 for every i. 

Now assume that there is no long root in a. We may assume that there 
is a root of the form e n _i — e n or e„_i + e n . 

In the first case, consider the basis a' of C n _i consisting of the roots 
{ei±ej £ a\l < i < j < n-l} and of the roots {e^ ± e n _i | ej ± e n G cr}. It 
is easy to see that the elements (<7i,<72, • • • ,<?n-i) are in T(a'). Indeed, for 
example, gig^-i = 1 if ± e n _i G a', as ft^ 1 = 1 and g n = g n ^ X - Thus by 
the induction hypothesis we obtain gf = 1 for all i 7^ n. Since c/„ = <? n _i, 
we also obtain gf l = l. 

The second case is similar. □ 



Let v = J27=i v i e i £ V be & vector with real coordinates and h = 
Y^=i hi e i £ V & vector with integral coordinates and such that Y17=i ^ 
is even. We will compute the normalized volume of Hc„(v) and the num- 
ber of integral points in Hc n (h) using Theorem 3.3. We will use the JK 
residue with respect to the measure dh associated to the basis ei, e2, . . . , e n . 
However, the normalized volume vo\z,c„(h) is computed for the measure 
determined by the lattice spanned by C n which is of index 2 in ®f =1 Zei. 

We introduce the function Jc{v) on U defined by: 

En 
i= l UiVi 

J c (v)(u) 



Ui=l 2u i ni<»<j<n(«i ~ U i) ]ll<i<j<n(^ + u i) ' 



For g = (51,52, . . . , g n ) G {±l} n the Kostant fraction (3) is the function 
on U defined by: 
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11=1 9 e^<=i 11 
Fc(g,h)(u) = — — '- 1 — u ._ u .)s 

1 

Theorem 9.2. Lei c be a chamber ofC(C n ). 

• For any v £c, we have 

vol z , c >) =2 JK C (J C («)). 

• For any vector h € V^n c uret/i integral coordinates such that Y17=i ^ 
is even, the value of the partition function is given by: 

Nc n (h)= Yl JK c (F c (g,h)). 

In the second formula, there should be a multiplication by a factor 2 as 
the volume of the fundamental domain of the lattice spanned by C n is 2. 
However, we should sum only on T = T/{±1}. Thus the two factors of 2 
compensate each other. In fact, we will indeed sum over T represented as 
{±l} n_1 x {1} and multiply the result by the constant 2. 

As in the case of B n , we will use the change of variable 1 + Zi = e Ui to 
compute more easily the formula for Nc n (h). As explained in the case of B n 
we need to use the integral formulation of the Jeffrey-Kirwan residue. Thus 
define 



c[9, )[ ) nr=i((i+^) 2 -i)xni< i <,<n(i+^-^(i+^)) 



i 



rii<i<j<n( 1 + z i)( l + z i) ~ 9i9j ' 

Performing the change of variables e Ui = 1 + Z{ on the function Fc(g, h)(u) 
and computing the Jacobian, Theorem 3.3 becomes: 

Theorem 9.3. Let c be a chamber ofC(C n ). 

• For any v 6 c, we have 

vol Z) c» = 2 JK C (Jc(«)). 

• For any vector h £ Vz H c with integral coordinates hi with Ya=i hi 
even, the value of the partition function is given by: 

NcM= E T^Tr^ L/c(g,h)(z)dz. 
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Similarly we will sum over T represented as {±l} n_1 x {1} and multiply 
the result by the constant 2. 

The cycle H(c) associated to a chamber c containing a regular element 
v = Y17=i Viei * s * ne same cycle that we computed in the preceding section 
for B n . Hence we can reuse most of procedures from the type B n . Paths 
are the same, and the residue calculations are the same. More precisely, the 
only two changes are in the computation of the set G (procedure GC(n)) and 
in the computation of the Kostant function (procedure UCKostant). This 
terminates the case of C n . 

10. The type D n 

10.1. The formulae to be implemented. Consider a vector space V with 
basis ei, . . . , e n . We choose the standard Lebesgue measure dh. Let 

T^n = {e i -e j \l<i<j<n}U{e i + e j \l<i<j< n}. 

Then V n is a positive roots system of type D n , and generates V. The number 
of elements in T> n is N = n 2 — n. 

We remark that the lattice Vz generated by roots of V n is the same 
lattice as the one generated by the roots of C n . It is of index 2 in L := 
Zei © Ze2 © • • • © Ze n and consists of elements v = [v± , v 2 , . . . , v n ] with 
integral coordinates such that the sum Y17=i v i ^ s an even integer. The 
group T = U /Uz is thus the quotient of T = U /Zei © • • • © Ze n , obtained 
by identifying g and —g, that is T = T/{±1}. As in Section 9, we identify 
the torus f = [//(Zei © ■ ■ ■ © Ze n ) = R n /Z n to (S 1 ) 11 by 

( Ul ,u 2 ,...,u n )^ (e 2 ^« V .. )e W^ _ 
Consider the set F = {±l} n C (S 1 )' 1 . For a a basic subset of V n , define 
f{a) = 1 5 e T e («- 2 ^v^TG> = j for alla G a j 

Lemma 10.1. Zei a be a basic subset ofV n . Then T{a) is contained in F. 

Proof. Basic subsets of V n are basic subsets of C n so that we can choose the 
same set F = {±l} n . □ 

Let v = Y17=i v i e i £ ^ be a vector with real coordinates and h = 
Y17=i hi e i ^ V a vector with integral coordinates and such that Y17=i 
is even. We will compute the normalized volume of IIx> n (v) and the number 
of integral points in U.x> n (h) using Theorem 3.3. 

Thus we introduce the function Jd{v) on U defined by: 

J D {v)(u) = 



Ili<i<j<„(«i - u j) Ili<i<j<„(«i + u j) ' 
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For g = (<?i,<72, ■ ■ ■ ,9n) £ {±l} n the Kostant fraction (3) is the function 
on U defined by: 



F D (9,h)(u) 



ni<i< j <„(i-0rV- (tti - w ' ) ) 
i 



X 



ni< i<j <n(l-5r 1 57 1 e- (ui+ ^ ) )' 



We have then 



Theorem 10.2. Lei c be a chamber ofC{V n ). 

• For any v £c, we have 

vohp n (v)=2JK c (J D (v)). 

• For any vector h € V^flc wi/i integral coordinates such that Y17=i hi 
is even, the value of the partition function is given by: 

Nv n (h)= £ JK C (F D ( 5 ,^)). 

ge{d=l}" 

We use the change of variable 1 + Zi = e Ui to compute more easily the 
formula for Nv n (h) and thus introduce integration over a cycle. Thus define 

TT n fl -I- 7^^+271-1-2 x rrn hi 



Ill<i<j<n( 1 + Z i ~ 9i9j(^ + Zj)) 
1 



After performing the change of variables e Ui = 1 + on the function 
Fu(g,h)(u) and after computing the Jacobian, Theorem 3.3 becomes: 

Theorem 10.3. Let c be a chamber ofC(T> n ). 

• For any v £ c, we have 

vol z ,©» =2 JK c (J D (u)). 

• For any vector h £ H c mi/i integral coordinates hi such that 
J2i=i hi i s even ; ^e value of the partition function is given by: 

N(V n ,h)= £ —^=- f T D (g,h)(z)dz. 

As for types A and B, in order to implement these formulae we first 
have to describe the set V(v,T> n ) (Section 10.2). We finish to explain the 
implementation of case D in Section 10.3, using the fact that types B and 
D are similar. 
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10.2. The search for maximal proper nested sets. A height function 
is 

n 

ht(v) = y](re — i)vi 
i=i 

which takes value 1 on all simple roots. We will deform it later on in order 
to have a function taking different values on roots. 

We now proceed to describe hyperplanes for D n . If P = [P + ,P~] are 
two disjoints subsets of {1, 2, . . . , n}, we denote by (up, v) the linear form 
^2ieP+ v i ~ J2jeP- v r Consider the hyperplane in V defined by 

H P = {v£ V,(u P ,v) = 0} 

and remark that it is equal to the hyperplane determined by the reverse list 
[P~,P + ]. Thus to each set P = {P + , P~} of two disjoint sets P + , P~ such 
that at least one is non empty, is associated a hyperplane Hp. 

We denote by Z the complement of P + U P~ in {1, 2, . . . , n} and by T>{Z) 
the subset of V n defined by 

V(Z) = {e i ±e j \l<i<j<n; i,j £ Z}. 

This is the positive roots system of type D\ z \, with the positivity induced 
by the lexicographic order. 

Let K.(P + , P~) be the subset of V n defined by 

{ei - ej | 1 < i < j < n; i,j £ P + } 
U { ei + e k \i £ P + ,k £ P-} 
U {e k - en | 1 < k < i < n; k,£ £ P~}. 

As we observed in Section 8.2 for B n , by defining /j = if i £ P + and 
fk = — e\p\-k+i if € P~ , the set K,(P + , P~) is a positive roots system of 
type A|p+| + |p-|_ 1 . Here the positivity is induced by the lexicographic order 
on P + and the reverse lexicographic order on P~ . 

Observe also that Hp is the vector space spanned by JC(P + ,P~) U V{Z). 

Lemma 10.4. • The hyperplane Hp is a V n - admissible hyperplane. 

• The set T> n n Hp is the union ofD(Z) and K,(P + , P~). 

• Every V n - admissible hyperplane is of this form. 

Proof. The first two assumptions are easy to see. Now as T> n is contained 
in B n , a P n -admissible hyperplane is B n admissible, so is of this form. □ 

10.3. The procedure MNS_KostantD. Most of procedures from type B n 
are kept unchanged. More precisely, the iterated residue calculation, the 
estimate of the order of poles and the global procedures coordinating com- 
putations are exactly the same as for type B n . 

The only serious adaptations to the case of D n appears in the procedure 
CheckDvector(n, v). In fact now we check that 
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for 1 < i < n — 1 , 
and is even, 
and is even. 

Other modifications are in procedures that are parent of CheckDvector. 
For example the procedure CheckDwall works exactly as CheckBwall, but 
now calls CheckDvector instead of CheckBvector. See Section 8.2. 

11. Performance of the programs 

In this Section, we describe several tests of our programs implementing 
the above MNS algorithms for types A n , B n , C n , D n . The algorithm im- 
plementation is made with Maple. Our programs are freely available at 
www . math . polytechnique . f r/~vergne/work/IntegralPoints . html. We 
compare our results with the ones obtained by two previous algorithms: 

• The Sp (for special permutations) algorithm by Baldoni-DeLoera- 
Vergne [2], only for A n . 

• The implementation LattE of Barvinok's algorithm [17], for every 
classical algebra; 

These two methods also helped us to test our algorithms on various ex- 
amples. 

Note that for our programs most of computation time is spent while 
computing iterated residues. Indeed MNS computation is fast and efficient. 
Note also that most of memory used by our programs serves to store all 
fractions that occur in the iterated residue process. The number of MNSs 
has a great influence on computation time, since we sum over all MNSs. 
In any case it seems that the deeper a vector is in the cone generated by 
positive roots, the higher the number of MNSs is. This is morally bound to 
the fact that there are more simplicial cones that might contain the vector. 
In Figure 10, we attach to every chamber c for B% the number of MNSs 
associated to any vector t)£c. 

Recall that the Sp algorithm relies on sums over a set Sp(a) of special 
permutations for a vector a. The main advantage of our algorithms is that 
we compute fewer iterated residues. In fact the number of MNS seems to 
be smaller than the number of special permutations that occur, for a given 
generic example. But, examples at the end of Table 11 show that a number 
of MNSs considerably smaller than those of Sp's doesn't lead to a better 
performance in time computation, even in the extreme case of just one MNS. 
Indeed this one residue computation can be very time consuming due to the 
substitutions z^ = Zj, which takes more time that the substitutions Zi = 
used in the Sp(a) algorithm. In the near future, we will improve this minor 
point. The MNSs method should be better and is better in general. 



vi + --- + Vi > 

vi H h v n -i + v n > 

vi -\ h v n -i -v n >0 
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e 2 -e 3 




e i-e 2 e 3 



Figure 10. Number of MNS containing any vector in a 
given chamber for B% 

During comparative tests, we figured out that one example in [2] has 
not been correctly copied from draft. More precisely in their Table 2 for 
complete graph K n , for the vector 

a = (82275, 33212, 91868, -57457, 47254, -64616, 94854, -227390) 

in the root lattice for A-?, the correct Kostant partition number is the 103- 
digits integer 

226040494681135377722281761934040091356424181 

242669497614801846058092972975120580334961426497 

and not only the first line of 45 digits. The Kostant number and Ehrhart 
polynomials for this a were computed on a 1GHz computer in 2,14 s and 
18, 54 s respectively, using 26 special permutations. Now with our programs 
running on a 1, 13GHz computer these times drop to 1,38 s and 2,50 s re- 
spectively, using 14 maximal nested sets. Similarly for the biggest example 
examined in [2], that is for the vector 

a = (46398, 36794, 92409, -16156, 29524, -68385, 
93335, 50738, 75167, -54015, -285809) 

in the root lattice for Aiq, the 189-digits answer was obtained in 2193 s using 
322 special permutations, whereas now we get the same result in 308 s using 
109 maximal nested sets. 

Table 11 contains respective performances for A n of LattE, Sp algorithms 
and our programs, a part the last four examples that compare only the last 
program with ours. Tables 12-14 contain respective performances for B n , 
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C n and D n of LattE and our programs. We also indicated the number of 
special permutations (Sp) and maximal nested sets (MNS). 

Tests were performed on Pentium IV 1,13GHz computers with 1500 or 
2000 mega-octets (Mo) of RAM memory. We stopped several computa- 
tions with LattE when we figured out that they would overcome computers' 
memory or take too much time with respects to the other algorithms; in this 
case we indicate the time spent and the number of mega-octets used by the 
computer. 
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60,4 s 
898 Sp 
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Root lattice element 


LattE 


MNS 


(1388,4024,3826) 


0,8s 


< 0, Is 
3 MNS 


(2691,5998,-6067,6184) 


2,6s 


0,1s 
1 MNS 


(1585, 7818, -2542, -2803, 2715) 
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3,0s 
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0,9s 
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(1082,947, 27, 42) 


22,9 s 


1,2s 
15 MNS 


(1047, 974, 20,44,-35) 


1939,9 s 


21,7s 
51 MNS 


(1015, 1082, -37, -21, -28, 14) 


> 7000 s 
> 1500 Mo 


378,0 s 
26 MNS 


Figure 12. Computation time for LattE and our programs, 
for B n 


Root lattice element 


LattE 


MNS 


(1388,4024,7652) 


0,8 s 


< 0, 1 s 
1 MNS 
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2,8s 
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0,8 s 


0,1s 
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12,2 s 


0,5s 
4 MNS 
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195,4 s 


2,8s 
6 MNS 


(1024,6,60, -6,-42,52) 


> 10800s 

> 2000 Mo 


1292,4s 
42 MNS 



Figure 13. Computation time for LattE and our programs, 
for C n 
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Root lattice element 


LattE 


MNS 


(8608, -305, 183) 


0,3 s 


< 0, 1 s 
1 MNS 


(32,5914,6166,-5360) 


1,5s 


< 0, Is 
1 MNS 


(1646, 3916, -3330, 6372, 7452) 


18,0 s 


0,5 s 
2 MNS 


(8127, 601, -2870, -2908, 10823, 3639) 


313,5s 


3,1s 
2 MNS 


(1009,1106,-9) 


0,2 s 


< 0, Is 
1 MNS 


(1074,959,64,77) 


3,0 s 


0,3 s 
6 MNS 


(1100,973,2,-1,-60) 


100,2 s 


3,1s 
18 MNS 


(1096,965,-54, 68,-34,-1) 


7076,7 s 


763,3 s 
47 MNS 



Figure 14. Computation time for LattE and our programs, 
for D n 



